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Abstract. Few-body physics has played a prominent role in atomic, molecular 
and nuclear physics since the early days of quantum mechanics. It is now 
possible — thanks to tremendous progress in cooling, trapping, and manipulating 
ultracold samples — to experimentally study few-body phenomena in trapped 
atomic and molecular systems with unprecedented control. This review 
summarizes recent studies of few-body phenomena in trapped atomic and 
molecular gases, with an emphasis on small trapped systems. We start by 
introducing the free-space scattering properties and then investigate what happens 
when two particles, bosons or fermions, are placed in an external confinement. 
Next, various three-body systems are treated analytically in limiting cases. Our 
current understanding of larger two-component Fermi systems and Bose systems is 
reviewed, and connections with the corresponding bulk systems are established. 
Lastly, future prospects and challenges are discussed. Throughout this review, 
commonalities with other systems such as nuclei or quantum dots are highlighted. 
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1. Introduction 

Although few-body physics plays an important role in all or nearly all branches of 
physics, few-body phenomena are often only covered briefly, if at all, in a typical 
undergraduate or graduate core curriculum. Only selected few-body problems such as 
the celestial three-particle Kepler problem and the quantum mechanical treatment of 
the helium atom, which becomes a three-body problem if the nucleus is treated as a 
point-particle, can be found in modern textbooks. One of the reasons is that few-body 
problems do not, in general, reduce to an effective one-body problem. Instead, the 
dynamics typically depend in non-trivial ways on all or nearly all degrees of freedom, 
making the analytical and numerical treatment of few-body systems in many cases 
challenging. This review discusses few-body aspects of trapped atomic and molecular 
Bose and Fermi gases, which are ideal model systems with which to study few-body 
physics. These systems exhibit a variety of few-body phenomena that are amenable 
to analytical treatments and that can, in certain cases, even be interpreted within an 
effective one-body framework despite the presence of strong interactions and multiple 
degrees of freedom. 

Degenerate Bose and Fermi gases are nowadays being produced routinely in 
laboratories all over the world. Most experimental set-ups work with alkali atoms 
(Li [H H], Na [3], K [1 [5], Rb [i [7], Cs 0), which have one outer valence electron. 
At ultracold temperatures, i.e., in the regime where the de Broglie wave length AdB 
is larger than the range tq of the two-body interaction potential, the interactions 
of the atomic gas are dominated by isotropic s-wave collisions. The net effect of 
an s-wave collision can be parametrized through the s-wave scattering length ao(0), 
which, in many instances, is the only microscopic parameter that enters into the 
description of larger systems [9l 110) . A well-known example of this type of description 
is the mean- field Gross-Pitaevskii equation [11] [T^ [T31 [HI [H], which depends on 
the product of the number of particles N (or more precisely, N — \ [16]) and the 
s-wave scattering length ao(0). The Gross-Pitaevskii equation has been proven to 
predict accurately a variety of aspects of dilute bosonic gases, including a mechanical 
instability for negative ao(0). Besides alkalis, atoms with a more complicated valence 
electron structure such as Ca [T^, He* [Tl[in], Cr [10], Yb [22122] and Sr [1S1[15 have 
been condensed. The atom-atom collisions between two Cr atoms, e.g., are anisotropic 
and long-ranged [25l [26] . These characteristics give rise to exciting beyond s-wave 
physics that depends on the strength of the dipole moment in addition to the s-wave 
scattering length [25l [26] . 

This review summarizes our current understanding of trapped few-atom systems 
that interact through isotropic short-range potentials. The quantized energy spectrum 
of the trapped system is the result of the interplay of the particle-particle interactions 
and the external trapping geometry. In fact, the external confinement can be viewed 
as quantizing the scattering continuum that is characteristic for the free-space system. 
Experimentally, trapped few-particle systems can be realized by loading so-called 
microtraps with just a few particles |27| . Alternatively, arrays of few-particle systems 
can be realized by loading a deep optical lattice, i.e., an optical lattice where tunneling 
between different lattice sites is strongly suppressed and where nearest-neighbor 
interactions can be neglected, with multiple atoms per lattice site [28l [29l [30] . 

A distinct characteristic of cold atom systems is that the atom-atom interactions 
can be tuned through the application of an external magnetic field in the vicinity of a 
Fano-Feshbach resonance [31] . A Fano-Feshbach resonance occurs when the two-body 
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energy of a closed channel molecule coincides with the scattering threshold of the 
open channel. Fano-Feshbach resonances in ultracold atomic gases were first observed 
experimentally in 1998 [311 [MUM]- Their existence paves the way to perform detailed 
comparisons between theory and experiment as a function of the interaction strength. 

Trapped few-atom systems — together with quantum dots, nuclei, atomic and 
molecular clusters, and certain nano-particles — fall into the category of mesoscopic 
systems. Quite generally, the study of mesoscopic systems provides a bridge between 
the microscopic and macroscopic worlds. This review discusses several links between 
small and large systems, and illustrates how the treatment of mesoscopic systems can 
provide an interpretation of the physics of large systems that is complementary to 
that derived within effective many-body frameworks. In this context, two-component 
Fermi gases with infinitely large s-wave scattering length provide a particularly rich 
example [351 133 133 [3S1 133- In this system, the underlying two-body interactions 
do not define a meaningful length scale, thereby making the system scale invariant. 
This scale invariance implies that the few-body system and the many-body system 
share the same or related universal features [40l |41] |42l |43l |44] . Moreover, the scale 
invariance implies that atomic two-component Fermi gases share certain characteristics 
with nuclear and neutron matter, despite the fact that the typical energy and length 
scales of these systems are vastly different. 

A particularly interesting few-body phenomenon, which is discussed in 
Subsecs. 13.31 14.31 and 15.31 is the Efimov effect. In the early 1970s, Vitaly 
Efimov [45l |46l considered a two-body system in free space with infinitely large s- 
wave scattering length that supports exactly one zero-energy bound state and no other 
weakly-bound states. Efimov found the surprising and, at first sight, counterintuitive 
result that the corresponding three-body system, obtained by adding a third particle 
and assuming pairwise interactions, supports infinitely many weakly-bound three-body 
states. Efimov's prediction, known as the three-body Efimov effect, stimulated much 
theoretical and experimental efforts [48l |49] , and continues to be one of the flagships 
of few-body physics today. Signatures of the three-body Efimov effect and associated 
Efimov resonances have been searched for in the helium trimer |50| and in nuclear 
systems [5T], among others. However, it was not till recently that the three-body 
Efimov effect has been confirmed conclusively experimentally via loss measurements 
of cold trapped atomic gases [52l [53l HH ESI EH E3 EH ES] ■ This review summarizes 
the key features of the Efimov effect, or more generally of Efimov physics, and its 
modifications when the three-body system is placed in an external trap. Extensions 
of the Efimov scenario to larger systems are also considered. 

The present review article aims at providing an introduction to the field of trapped 
few-atom and few-molecule systems. Readers who would like to place this topic into a 
broader context are directed to a number of other review articles. Experimental and 
theoretical aspects of cold and ultracold collisions have been reviewed by Weiner et 
al. in 1999 [50] ■ The topic of Fano-Feshbach resonances has been reviewed by Chin 
et al. in 2010 [31] (see also Refs. [61] [62]). The Efimov effect in free-space has been 
reviewed by Braaten and Hammer in 2006 and 2007 [48l|49]. A related "Trend" was 
authored by Ferlaino and Grimm [63], and a Physics Today article by Greene [6^ in 
2010. The broad subjects of Bose and Fermi gases are covered in reviews by Dalfovo 
et al. from 1998 [TS] and Giorgini et al. from 2008 [31], respectively. Many-body 
aspects of cold atom gases have been reviewed by Bloch et al. in 2008 [52]. Various 
reviews have focused on specific few-body systems: Jensen et al. discuss the physics 
of halo systems [66] , Carlson and Schiavilla that of few-nucleon systems [67] , Reiman 
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et al. that of vortices in few-particle droplets [5^ , and Nielsen et al. that of three- 
body systems with short-range interactions [51] . Lastly, the treatment of few-particle 
systems within the hyperspherical framework, which is utilized at various places in 
the present work, has been reviewed by Lin in 1995 |70| and Rittenhouse et al. in 

2011 [H]. 

The remainder of this article is structured as follows. Section [2] discusses 
the underlying two-body physics, both without and with confinement, while Sec. |3] 
considers the physics of three interacting particles under external spherically 
symmetric confinement. Section 2] discusses various aspects of trapped Fermi gases, 
focusing primarily on s-wave interacting two-component Fermi gases. Section [5] 
summarizes the key characteristics of s-wave interacting Bose gases under external 
confinement. Lastly, Sec. [H] comments on a number of present and future frontiers. 

2. Underlying two-body physics 

2.1. Free-space two-particle scattering in three-dimensional space 

This subsection considers the two-particle system in free space while the next 
subsection treats the trapped two-particle system. For two neutral atoms, the 
two-body potential at large interparticle distances r is dominated by the van der 
Waals tail that falls of as —Ce/r^ (Cg denotes the van der Waals coefficient). The 
characteristic length or range /char of the two-body interaction potential is defined 
as the distance at which the kinetic and potential energies are approximately equal, 
/char — (S/iCg/fi.^)^/^ 60 , where /i denotes the reduced mass of the two interacting 
particles. For alkali atoms, /char is of the order of lOOabohr, where flbohr denotes the 
Bohr radius, Cbohr — 5.292 x 10~^^m. 

The potentials that describe diatomic alkali systems support many two-body 
bound states in free space. For example, the number of vibrational bound states in the 
rotational ground state is respectively 66 and 16 for the electronic ground and excited 
state X^E+ and a'^I]+ potential curves of ^■^Na2 [75]. Throughout this review, we 
are concerned with the low-temperature regime, i.e., we are interested in physics that 
occurs at energy scales that are much smaller than i^char, where £^char = ^^/(2A*^char)' 
or equivalently, length scales that are much larger than /char- This includes at most 
a few weakly-bound two-body states as well as the energetically low-lying scattering 
continuum (or "quantized scattering continuum") but excludes deeply bound two- 
body states and the high-energy scattering continuum. Generally speaking, the low- 
energy physics is independent of the details of the underlying two-body potential and 
is determined by a few parameters such as the two-body scattering length and the 
effective range. Thus, low-energy phenomena can be described by replacing the true 
or full interaction potential by a model potential that is characterized by the same 
parameters as the true potential (e.g., the same scattering length and effective range), 
or by effective theories such as an effective range expansion approach or an effective 
field theory approach. Effective range and effective field theory approaches have been 
applied to a variety of systems, ranging from the scattering of atoms at ultracold 
temperatures [15] to the low-energy scattering of two nucleons in free-space [75] and on 
a lattice [TH [751 [75] to the description of particle physics phenomena such as narrow 
charmonium resonances [77]. The effective pseudopotential description, frequently 
employed in this review, goes back to Fermi, who treated the scattering between slow 
neutrons and hydrogen atoms |78j . 
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In the following, we assume that the two particles of mass mi and TO2 (with 
position vectors ri and interact through a spherically symmetric two-body 
potential 14b (^) 7 where r denotes the distance coordinate between the two particles, 
r = |ri — r2|. Separating off the center-mass degrees of freedom, the relative scattering 
wave function i/'fe(r) can be decomposed into partial waves, i^kir) = X)/ Rik{f')Pi{'d), 
where denotes the polar angle and Pi the Legendre polynomial of degree I. 

Asymptotically, i.e., in the regime where the interaction potential can be 
neglected, the radial wave functions are given by [TS] 

R°k\r) = Ni{k) {ji{kr) ~ tB.n[5i{k)]ni{kr)} , (1) 



where k is related to the scattering energy E through k = \/2fiE / h and where ji and n; 
denote the spherical Bessel function and the Neumann function, respectively. In the 
small kr limit, ji{kr) and ni{kr) behave as (fcr)' and (fcr)~('+^\ respectively |79) . 
Correspondingly, ji and ni are referred to as regular and irregular solutions. In 
Eq. ([1]), Ni{k) denotes a normalization constant. The scattering phase shifts Si{k) 
are determined by matching the asymptotic solution i?™*(r), Eq. (HI), to the inner 
solution R"^{r), which depends on the details of the potential Vth{r). The phase shifts 
accumulate in the inner region and quantify the admixture of the irregular solutions 
ni{kr). A vanishing phase shift 5i indicates an effectively non- interacting system, 
i.e., the asymptotic solution of the 1*^ partial wave channel is identical to that of 
the corresponding non-interacting system (differences can exist, though, in the inner 
region). A small positive phase shift indicates that the asymptotic solution of the 
Z*'* partial wave channel is shifted to smaller kr ("pulled inward") compared to that 
of the non-interacting system, thus indicating an effectively attractive interaction. A 
small negative phase shift, in contrast, indicates that the asymptotic solution of the 
l^^ partial wave channel is shifted to larger kr ( "pushed outward" ) compared to that 
of the non-interacting system, thus indicating an effectively repulsive interaction. The 
phase shifts 6i{k) can be used to define the energy-dependent generalized scattering 
lengths af'^^{k) [79] (note that the opposite sign convention is sometimes employed 
in nuclear and particle physics), 

«rW = -^^- (2) 

For short-range potentials, i.e., for potentials that fall off faster than l/r^'+s 
asymptotically, the energy- dependent scattering lengths approach a constant in the 
zero-range limit [HD], a^'^^(O) = Vmvk^Q af''^^ {k) (for long-range potentials, see 
Subsec. 12. 3p . For the lowest partial wave, i.e., for I — 0, ao(0) is the usual s-wave 
scattering length, which dominates, away from higher partial wave resonances, the 
low-energy cross sections of bosonic gases, two-component Fermi gases and Bose- Fermi 
mixtures. In the zero-energy limit, the scaled outside wave function ri?QQ*(r) becomes 
A'^o(O) [r — ao(0)]. The leading-order energy-dependence of the s-wave scattering 
length can be written in terms of the effective range feff [79], 

— TTT ~ - ^^offfc^; (3) 

ao(/c) ao(0) 2 

the validity of Eq. ([3]) requires that the potential falls off faster than 
asymptotically. For identical fermions, the Pauli exclusion principle forbids s-wave 
scattering, making the p-wave channel with generalized scattering length or scattering 
volume a\{k) dominant. 

Figure [Tj illustrates the behavior of the scattering lengths and wave functions for 
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Figure 1. (Color online) Illustration of two-body scattering for an attractive 
square well model potential with depth Vb and range tq. (a) Solid and dotted 
lines show the zero-energy scattering lengths ao(0) and af (0) as a function of 
'^in^'Oi where k?^ = 2^Vo/h^. Inset: Solid lines show the energy of the most 
weakly-bound s-wave state for the square well potential as a function of k\^rQ. 
For comparison, dotted lines show Eq. jTj with ao(fc) approximated by ao(0) of 
the square well potential, (b) Solid, dotted and dashed lines show the scaled 
s-wave scattering wave function rRQo{r) as a function of r/ro for [ao(0)]~^ = 
(one zero-energy s-wave bound state), ao(0) = — ro/2 [no s-wave bound state; 
ao(0) is marked by a circle] and ao(0) = ro/2 [one s-wave bound state; ao(0) is 
marked by a square], respectively. 



an attractive square well potential with depth Vb and range ro- This or similarly 
simplistic model potentials suffice to describe generic low-energy features of cold atom 
systems. Solid and dotted lines in Fig.[IJa) show the scattering lengths ao(0) and a\{Q) 
as a function of the dimensionless "strength parameter" kiy^r^, where fcm = VS/^Vb/fi. 
The scattering lengths ao(0) and 0^(0) vanish for vanishing Vb, and decrease with 
increasing kinTo- The s- and p-wave scattering lengths diverge when ki^rQ = 7r/2 and 
kinfo — TT, respectively. As /cinT-Q increases further, the scattering lengths decrease 
from large positive to large negative values. When the scattering lengths ao(0) and 
(0) diverge, the square well potential supports a new s-wave or a new p-wave bound 
state, respectively. Dashed, dotted and solid lines in Fig. [TJb) show the scaled s- 
wave scattering wave function rRooir) for finite positive, finite negative and infinitely 
large ao(0), respectively. The linear behavior of r_Roo('') for r > ro is clearly visible. 
For finite ao(0), the scattering length can be "read off" by extrapolating the linear 
behavior to r < tq [see symbols in Fig. [Ijb)]. By construction, the wave functions 
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for the same zero-energy scattering length but different depth or number of bound 
states (not shown) agree in the large r limit. Differences, however, exist in the small 
r region. These short-range features are typically not probed at low temperatures, 
making the scattering length the only relevant quantity. This can be intuitively 
understood by realizing that the de Broglie wave length AdB increases with decreasing 
temperature, thereby setting a resolution limit, i.e., the collision energy (which is 
in typical experiments set by a combination of the temperature and the particle 
statistics) is too small to probe features at small length scales. Correspondingly, 
the low temperature regime is known as the long wave length limit. 

Mathematically, it can be advantageous to replace the true interaction potential 
or shape-dependent model potential by a zero-range pseudopotential that acts only 
when the two particles sit on top of each other [7S1 IHOl HI]- Using the quantum 
defect language (see, e.g., Refs. [551 [531 |H3), the pseudopotential can be thought of 
as enforcing a particular value of the logarithmic derivative at r = 0. In three spatial 
dimensions, the regularized s-wave Fermi-Huang pseudopotential V^"^ reads [81] 

K)''M=5o(fc)<5(r)|:r, (4) 

where the interaction strength go{k) is chosen such that V^^ reproduces the energy- 
dependent s-wave scattering length ao{k) of the true interaction potential, 

9o[k) = . (5) 

The regularization operator {d/dr)r in Eq. Q ensures that the pseudopotential is 
well behaved as r goes to [81]. The pseudopotential imposes a boundary condition 
on the relative wave function ipkir) at r = [85» S&i , 



^ (6) 



ao(fc) ■ 

^0 

This so-called Bethe-Peierls boundary condition serves as an alternative parametriza- 
tion of VQ^{r). The pseudopotential Vq^ supports no two-body bound state for neg- 
ative s-wave scattering length ao(fc) and a single two-body s-wave bound state with 
binding energy inbound, 

" 2/i[ao(fc)]2' 

for positive ao(fc). Equation ([7|) is an implicit equation for the bound state energy, 
which requires knowing ao(fc) for imaginary k. The size of the bound state of Vq^ 
is determined by ao(fc), i.e., ^ (r^) = ao{k)/^/2. We note that Eq. ([7]) also describes 
the s-wave bound states of finite-range potentials. This can be seen by analytically 
continuing Eq. ([1]) to negative energies or by analyzing the poles of the energy- 
dependent S-matrix [79]. In the low-energy limit, Eq. ([7]) provides a good description 
of the bound state energies if ao(A:) is approximated by ao(0) (see also below). In 
fact, it should be noted that the original formulation of the pseudopotential (or 
corresponding boundary condition) employed the zero-energy scattering length oq (0) 
instead of ao(fc). In the cold atom context, the energy-dependent scattering length 
ao{k) was first introduced in Refs. [571155] . 

It is important to note that V^''^ only acts on the lowest partial wave. The 
partial waves Rfl^^ with / > 1 have vanishing amplitude at r = 0, and are thus not 
affected by Vq^. This implies that Eq. ([6]) remains valid if ipk is replaced by i?Q^* and 
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that Rf^*'{r) = Ni{k)ji{kr) for / > 1. Pseudopotentials VJ^^ for higher partial wave 
scattering have been developed in the literature [HIllHniliniinillMIMllSlinS- Similar 
to the s-wave pseudopotential, the higher partial wave pseudopotentials, applicable to 
spherically symmetric systems, can be parameterized in terms of a boundary condition 

at r = 0, 

q2!+l 



r(z + i/2) 
V^Tii + 1) 



I^TTT {r'+^ J ^fc(f)Pi (cos cos ^) 
J ipi^(^f) Pi (^cos^)d cos 



Solid lines in the inset of Fig. [IJa) show the energy of the most weakly-bound 
s-wave state of the square well potential as a function of hnro- As discussed in 
the context of the main figure, a new s-wave bound state is being supported when 
hnro = mr/2, where n = 1, 3, • • •. For comparison, dotted lines show the bound state 
energy -Ebound, Eq. ([7]), with ao(fc) approximated by the zero-energy scattering length 
ao(0) of the square well potential. For \E\/Echai ^ 1, where E'char = f^/i'^t^rl) (i.e., 
in the low-energy regime), the agreement between the solid and dotted lines is good. 
However, as l-Ej/iJchar increases, deviations between the two data sets are visible. This 
can be attributed to the fact that the energy-dependence of ao{E) becomes stronger 
as increases. 

Typical interactions between alkali atoms involve two electronic energy curves, 
the singlet and the triplet Born-Oppenheimer (BO) potential curves labeled X 
and a '^Sj, respectively, which are coupled through a hyperfine Hamiltonian. The 
coupling is responsible for the multi-channel nature of the scattering process. Multi- 
channel scattering is reviewed in detail in Refs. [601 US- For our purposes, it is 
important that the coupling between the channels can give rise to Fano-Feshbach 
resonances |31j . A Fano-Feshbach resonance in the Ith partial wave arises when the 
scattering energy coincides with the energy of a bound state in the Ith partial wave. 
Since the "detuning" between the open and closed channels in the large r regime is 
determined by the difference in magnetic moments, the thresholds of the two channels 
can be varied through the application of an external magnetic field. If the topology of 
the potential energy surfaces is such that the variation of the external field allows for 
the bound state energy and the scattering energy to coincide, then the system exhibits 
a resonance. Varying the strength B of the external magnetic field in the vicinity of a 
Fano-Feshbach resonance allows the coupled-channel Z-wave scattering length af^'^^^k) 
to be tuned to essentially any value. 

Assuming no overlapping Fano-Feshbach resonances exist and higher partial 
wave contributions are negligible, the magnetic field-dependent coupled-channel s- 
wave scattering length ao{B) can be parameterized by three effective parameters, the 
background scattering length abg, the resonance width AB and the resonance position 

-Bres [22], 

ao(S)= abg (^1-^^41;-). (9) 

This effective description has been confirmed by coupled-channel calculations and 
experimental data; it can, for example, be derived using a simple coupled-channel 
zero-range or square well model [96[ 197] . The parameters A_B, B,-cs and Obg are 
extracted from experimental data or full coupled-channel calculations. The solid line 
in Fig. [2ja) shows the magnetic field dependence of the scattering length for the 
magnetic Fano-Feshbach resonance for two ^'^Rb atoms in the {F, Mp) = (2, —2) state 
near 155G (here, F denotes the hyperfine quantum number and Mp the corresponding 
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Figure 2. (Color online) (a) The solid line shows the field-dependent scattering 
length ao(-B), Eq. for the Fano-Feshbach resonance around 155G for two *^Rb 
atoms in the (_F, Mp) = (2, —2) state as a function of B\ the parameters adopted 
are A_B = lOG, -Bres = 155. 2G and a^g = — SSOa^ojir |97l . For comparison, 
the dotted line shows the effective scattering length for a "mock resonance" with 
AS = IG but the same abg and Bres as for the *^Rb resonance, (b) Solid and 
dotted lines show the s-wave energy spectra, calculated using Eqs. I I17I1 and | [9t . 
for two atoms under harmonic confinement with ajj^^^ = 5000ajjoijj, . 



projection quantum number). It is often useful to write the field-dependent scattering 
length in terms of the energy width T e |72| , 

where Te — 2kahgAB[dE^os{B)/dB] and where Ercs{B) denotes the field-dependent 
bound state energy of the coupled-channel system. Here, E and k denote, as above, 
the two-body scattering energy and corresponding wave vector, respectively. We 
note that the £^-term in the denominator of Eq. pUj) can often be neglected. In 
this case, Eq. (|T0)) reduces to Eq. ([9]). The quantity dErcs{B)/dB is known from 
experiment and/or theory. The energy width Te can be used to define broad and 
narrow resonances [Ml IMl ESI HOQl [Ml [l02]. If is larger (smaller) than the 
characteristic energy scale of the system, the resonance is called broad (narrow). 
The characteristic energy scale depends on the system under investigation. For the 
BCS-BEC crossover problem, for example, the characteristic energy scale in the low 
temperature regime is set by the Fermi energy when |ao(-B)| is small and ao{B) < 0, 
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and by the two-body binding energy when ao{B) is small and ao{B) > 99J. As a rule 
of thumb, broad resonances can be modeled quite accurately within a single-channel 
model; narrow resonances, in contrast, require a coupled-channel treatment. 

2.2. Two particles in a spherically symmetric harmonic trap 

This subsection considers two particles interacting through a spherically symmetric 
potential Vtb under external confinement, i.e., we consider what happens to the 
scattering continuum and the weakly-bound states when the asymptotic boundary 
conditions are modified by the external confinement. Throughout, we consider 
the regime where the confinement length is much larger than the range of the 
underlying particle-particle interaction potential. Experimentally, trapped systems 
can be realized by loading atoms into an optical lattice created by counterpropagating 
standing light waves. The atom-light interaction creates a periodic potential felt by 
the atoms. Near the bottom of each well, the confining potential is approximately 
harmonic. Furthermore, for atoms with short-range interactions, the spacing between 
the lattice sites is so large that "off-site" interactions can be neglected. Assuming 
that tunneling between neighboring lattice sites can also be neglected, the lattice sites 
are independent from each other and each lattice site hosts an independent few-body 
system. Experimentally, the number of atoms per lattice site can be controlled. Here, 
we treat the two-body system; larger systems are discussed in Sees. 

We assume that the external confinement is harmonic, i.e., that the i*^ particle 
feels the external potential miUj'^ff /2 (i = 1 or 2), where uj denotes the angular 
trapping frequency and r-i the position vector of the i*^ particle measured with respect 
to the trap center. For this system, the center of mass motion separates, i.e., -Etot = 
Ecm + E, where the center of mass energy is given by Ecm — (2(3cm + Lcm + 'i/2)huj 
with Qcm = 0, 1, • • • and icm = 0, 1, • • •, and E denotes the relative energy. Assuming 
a spherically symmetric interaction potential Vth{r), the relative Hamiltonian for the 
two-body system reads 

^^ = -^v|+iAic.V + ytbM, (11) 

where /x denotes, as above, the reduced two-body mass. The relative eigen functions 
ipiq{r) separate, i^iq{r) = r^^uiq{r)Yijni'&, 'fi). uiq{r) solves the radial Schrodinger 
equation 

Ulq{r) = ElqUlq{r), (12) 

where Eiq — E ^ {2q + / + 3/2)huj. The non-integer quantum number q labels the 
radial solutions for a given relative orbital angular momentum /. Equation (I12p can 
be solved numerically for essentially any Vth{r). Here, we seek analytical solutions for 
the l-wave zero-range pseudopotentials Vi^^. 

The linearly independent solutions fiq and giq of the radial Schrodinger equation 
with Vth = Vi^^ for r > are identical to those for an isotropic harmonic oscillator 
with mass ^. Introducing the harmonic oscillator length aho,;i, Oho.p — \/^7{p^-, and 
defining the dimensionless variable x, x = r/aho,p,, fiq and giq can be written as 

fiqix)=x'+'e^p(^) iFi(-<7,/ + 3/2,x2) (13) 



-h-" hH{l + 1) 1 



2/i dr'^ 



2/ir2 
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and 

gi^{x)^x~'e^-p{^^ iFi{-q- I- 1/2,-1 + l/2,x^), (14) 

where iFi denotes the confluent hypergeometric function of the first kind. We write 
uiq{x) = ci{q)fiq{x) + di{q)giq{x) and choose the ratio ci{q)/di{q) such that uiq{x) 
vanishes in the large x hmit, i.e., so that the exponentially increasing pieces of fiq{x) 
and giq{x) cancel as a; — ^ oo. Using the a; — >■ oo behavior of iFi [103] . we find 

r(/ + 3/2) rH + 1/2) 

T{-q) T{-q- 1/2- I) 

Using Eq. p5|) in uiq{x), we find 



uiqix) ^ Niiq) exp (^-^^ x'+'U{-q, I + 3/2, x^). 



(16) 



where U denotes the confluent hypergeometric function of the second kind and Ni (q) 
a normalization constant. Next, we ensure that r~^uiq obeys the boundary condition 
in the small x limit, i.e., we enforce Eq. The resulting quantization condition 
reads 



(■ 



2hu} 2 ' 4 



The transcendental equation, Eq. (IT7)) . determines the eigenspectrum for any I and 
af''~^^{k). The eigenenergies Eiq can be obtained straightforwardly using a simple 
rootfinding routine. For 1 = 0, this eigenequation was first derived by Busch et al. [104) . 
For / > 0, the solutions can be found in Refs. [HllIIl UlSllini] . The wave function 
uiq is nornializable for / = but not for I > 0. This issue has been discussed in detail 
in Refs. [95, 106, 107| and will not be elaborated on here. 

Figures 131^ a) andlljb) show the energy spectra for / = and / = 1, respectively, 
as a function of the inverse scattering lengths. Solid lines show the energy spectra 
obtained from Eq. (|17p using a^'^^(fc) = af^~^^{0). For comparison, dotted lines 
show the eigenspectra obtained for a square well potential with a fairly large range 
''0: ''0 = 0.1aho,/i. The s-wave spectrum for the square well potential is fairly well 
reproduced by the zero-range model for all ao(0); the deviations are largest for negative 
energies where the short-range physics plays a non- negligible role. For 1 = 1, the square 
well spectrum is qualitatively reproduced for E > 5huj/2 and large |aho,/^/ii(0)|^. In 
other parameter regimes, significant deviations are visible. A significantly improved 
description is obtained by treating the energy-dependence of the generalized scattering 
length explicitly. To this end, we solve Eq. (|17l) self-consistently f87l [88], i.e., we 
seek solutions for which the energy on the left hand side of Eq. (fT7|) is the same 
as that at which the generalized free-space scattering length af'^^(fc) for the square 
well potential with = O.laho,^ is being evaluated. The resulting eigenenergies [not 
shown in Figs.[31[a) and[31[b)] are in excellent agreement with the eigenenergies for the 
square- well potential for both I = and 1 = 1. In general, it is found that the use of 
the energy-dependent scattering length is crucial for the description of higher partial 
wave physics within the zero-range pseudopotential framework [911 192[ [51] . 

Figure [H shows the / = wave function itoq(r) for several ao(0) corresponding to 
the lowest energy branch. In the small |ao(0)| limit, ao(0) < 0, uoq{r) approaches the 
non-interacting wave function. In the small ao(0) limit, ao(0) > 0, in contrast, uoq{r) 
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Figure 3. (Color online) Illustration of the energy spectrum of the trapped two- 
particle system for (a) 1 = and (b) 1 = 1. Solid lines show the eigenenergies 
Eiq obtained by solving Eq. I I17II with a^'^^(fc) approximated by a^'^^(O). Dotted 
lines show the eigenenergies obtained for a square well potential with rg = O.laho.ji 

(0) [and such that the free-space 



,21-1-1 



and Vb adjusted to obtain the desired 
system supports no (one) Z-wave bound state for a^'^^(O) < (a^''''^(0) > 0)] 




Figure 4. (Color online) Dotted, solid and dashed lines show the relative I = 
wave function uoq{r), Eq. HIGH , corresponding to the lowest energy branch of the 
trapped two-body system with ao(0)/aijo,(j = —0.1, oo and 0.1, respectively. 
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approaches the wave function of a tightly bound molecule in free-space. At unitarity, 
the admixture of the irregular solution is maximal, leading to a wave function that has 
vanishing derivative at r = 0. Formally, this wave function is identical to the relative 
even-parity wave function of two non-interacting harmonically trapped particles in 
one dimension. 

We now discuss the dependence of the s-wave energy spectrum on the external 
magnetic field strength B. Solid and dotted lines in Fig. [Sfa) show the magnetic 
field dependent scattering length ao(i3), Eq. (|9]), for two resonances with the same 
abg and Bros but different widths AS, i.e., Ai? = lOG and IG, respectively. Solid 
and dotted lines in Fig. [Sfb) show the s-wave energy spectra obtained from Eq. (ITTl) 
corresponding to the resonances shown in Fig.ISfa). To this end, the scattering length 
ao(0) in Eq. pT|) is converted to B via Eq. Figure [UJb) shows that the energy 
spectrum for A_B — lOG is characterized by broader avoided crossings than that for 
AB = IG. The coupling between the resonance state (or the closed molecular channel) 
and the trap states increases with increasing /S.B. A Landau-Zener analysis [7^ 
shows that the energy width F^, Eq. p^ . determines the probability to make a 
transition from the trapped state to the molecular state. In particular, it is found 
that the transition probability decreases exponentially with increasing F^; for the 
same ramp speed. Correspondingly, the closed molecular channel plays a negligible 
role for resonances with relatively large Ai3 and a non-negligible role for those with 
relatively small AS. Thus, broad resonances can, generally, be described within a 
single channel framework while narrow resonances require a two-channel description 
(see also Subsec. 12. ip . 

2.3. Beyond spherical symmetry 

The previous subsection discussed the physics of two atoms interacting through 
a spherically symmetric potential under spherically symmetric confinement. This 
confining geometry is fairly easy to treat both numerically and analytically, and is 
directly relevant to on-going cold atom experiments. This subsection considers systems 
for which the spherical symmetry is broken either due to a non-spherically symmetric 
confinement, a non-spherically symmetric interaction potential or both. In all three 
cases, new two-body resonances arise that can be tuned by varying the confinement 
(thus motivating the term confinement-induced resonance [108[ 1109] ) or by changing 
the admixture of the non-spherical contribution to the interaction potential. These 
tunable two-body resonances pave the way for investigations of intriguing few- and 
many-body phenomena. 

Table [1] lists confinement geometries that have been considered in the literature. 
As an example, we consider two particles interacting through a spherically symmetric 
potential Vtb under cylindrically-symmetric confinement in more detail |94[ I108[ 
[TMl [TTnl [iTTl [TT21 [TT3l [TT4] . As in Subsec. [221 the center-of-mass degrees of 
freedom separate off and the confinement in the relative coordinate becomes T4rap = 
/x(A^aj^p^ +w^z^)/2, where z = zx — Z2 and p — {x\ — X2Y + {yi — 2/2)^- The aspect 
ratio A is defined in terms of the angular trapping frequencies cvp and ujz along the p- 
and z-directions, A — LUp/cuz. Intuitively, one expects that a tight confinement in one 
spatial dimension (A <C 1) leads to an effectively two-dimensional system in which the 
motion along the tight confinement direction, or high energy coordinate, is frozen out. 
Correspondingly, a tight confinement in two spatial dimensions (A ^ 1) is expected 
to lead to an effectively one-dimensional system. This argument can be formalized by 
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Table 1. Selected trapping geometries considered in the literature for two-particle 
systems. Other geometries such as double well geometries and optical lattice type 
potentials have also been considered. 



Trapping geometry 


References 


hardwall at r = R (three-dimensional) 
^liuj'^r^ (three-dimensional) 
^fiuj^p^ (two-dimensional) 
■^IMaPz^ (one-dimensional) 
|/i(a;2p2 j^uilz^) 

i(miajfrf -|-m2aj^rf) 

cubic box; periodic boundary condition 

box with lengths Lx, Ly, L^; periodic boundary condition 


m 

[841 Ion [921 [Toil [T05] 
[1151 

[Ml[T08l[T09l[TTol[TTTl[TT2l[TTl[TT4] 

[TT611TT711TT8) 

[TT911T2011T2T] 

EH 



defining efi^ective one- and two-dimensional interaction potentials. 

We exemplarily consider the A ^ 1 case with s-wave interactions [1081 1109) . 
To start with, we assume that the motion in the p- and z-directions separates and, 
further, that the motion along the p-coordinate is frozen in the harmonic oscillator 
ground state 0ho,o(p)- Modeling the two-body interaction potential through V^'^ , we 
define an effective one-dimensional potential 1/3^(2;) through 

/•oo 

V^i{z)^2n ct>LAPWo''{^<l>i^oAp)pdp- (18) 
Jo 

The integral evaluates to V^ff^(2;) = glgS{z) with g^^ — 2h'^ai^{Q) / {pLa\^ p), where 

<i\lo,p denotes the transverse confinement length, aho,p — \/h/ {p,tOp)- So far, we have 
assumed that the motion along the p-coordinate is frozen in the harmonic oscillator 
ground state 0ho,o(p)- This is, however, not the case since the pseudopotential Vq^ 
couples the p and z degrees of freedom. Thus, the cylindrically symmetric outside 
solution (valid for r > 0) needs to be projected onto the spherically symmetric inside 
solution (valid for r = 0), thereby leading to an admixture of excited state harmonic 
oscillator wave functions 0ho.j(/3); here, j collectively denotes the quantum numbers 
associated with the motion in the xj;-plane. For finite-range potentials, the outside and 
inside solutions can, using a frame transformation approach, be connected at a finite 
r value |112l . A full microscopic calculation for the zero-range potential gives |108| , 



ffeff — 2 



(19) 



^ho,p 

where C(l/2) —1.4604. For smaU |ao(0)|/aho,p, Eq. ((T^ reduces to the naive 
estimate discussed above. For large |ao(0)|/aho,p, however, the strong transverse 
confinement modifies the effective one-dimensional coupling constant. In particular, a 
diverging one-dimensional coupling constant, and thus a one-dimensional resonance, 
is found for ao(0)/aho,p = |C(1/2)|~^- Effectively one-dimensional many-body systems 
are, in many cases, well described by strictly one-dimensional Hamiltonian such as 
the Lieb-Liniger or Gaudin-Yang Hamiltonian |123[ I124[ 11251 1126) . For cold atom 
realizations of these Hamiltonian, the coupling constant is given by g^g, Eq. p^ . 
Since g],^ can be tuned by varying aho,p or ao(0), Eq. forms the basis for 

studying strongly correlated one-dimensional systems. One-dimensional systems have 
attracted a great deal of attention from mathematicians and physicists since they 
are more amenable to analytical treatments than two- or three-dimensional systems. 
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Strictly one-dimensional systems with contact interactions, e.g., are integrable. 
Deviations from a strictly one-dimensional geometry may lead to deviations from 
integrability. Thus, important questions regarding integrability, near-integrability 
and thermalization can be investigated by studying highly-elongated quasi-one- 
dimensional Bose and Fermi gases [123 [US [US ITBOlfm] . 

Similar analyses have been conducted for higher partial waves as well as for 
geometries with A ^ 1 [94 | 11101 1112j . In all cases, confinement induced resonances, 
i.e., diverging effective low-dimensional coupling constants, are found for finite 
three-dimensional scattering lengths. The study of small quasi-two-dimensional 
trapped Fermi gases is partially motivated by analogous studies of two-dimensional 
quantum dot systems |132[ 1133) . In the context of cold Fermi gases, the filling of 
energy levels has been analyzed using Hund's rules [134] . An alternative trapping 
geometry consists of a periodic optical lattice, which also gives rise to confinement 
induced resonances that can be tuned by varying the lattice height or other system 
parameters [HSl [^1 [TOTl ITBSj . 

The energy spectrum of two-particle systems in a cubic box with periodic 
boundary conditions has been analyzed extensively in the context of lattice QCD 
calculations [74j [75j 1139] . For this confining geometry, the center of mass of the 
system does not decouple and two-body interactions lead to a modification of the 
plane wave solutions applicable to the non- interacting system. In analogy to Eq. pT)) . 
the two-body energy spectrum for the box- confinement has been related to the phase 
shifts. The corresponding equation is termed Liischer's formula [75] and has found 
applications in QCD and other lattice simulation studies. An analysis of the energy 
spectrum determined from Liischer's formula shows that the relationship between the 
phase shifts and eigenenergies is notably different from that found for the spherically 
symmetric harmonic confining potential treated in Subsec. [2.2l This shows exemplarily 
that the energy spectrum depends in non-trivial ways on the geometry of the external 
confinement. 

Lastly, we comment briefly on two particles that interact through a non- 
spherically symmetric dipole-dipole potential Vdd under spherically symmetric 
confinement. For simplicity, we assume that the electric or magnetic point dipoles are 
aligned along the z-axis, Vdd{r) = d^(l — 3cos^ '&)/r^, where d denotes the strength of 
the dipole moment and 'd the angle between the z-axis and r. For r < tq, we assume 
a spherically symmetric potential Vsr{r). The dipole-dipole potential introduces a 
new length scale, the dipole length d*, rf* — 2fid^/h'^. Since the angle-dependence 
of the dipole-dipole potential leads to the coupling of different partial waves, the 
scattering framework discussed in Subsec. 12.11 requires modifications. One of the 
consequences is that one now has to consider a scattering matrix with elements 
am, where aiq{k) = - ta,n[Si'i{k)]/k [60 } 11401 fT4T| 1142] . Compared to short-range 
potentials [see Eq. ([2|)], the dipole-dipole interaction potential leads to a modified 
threshold behavior, i.e., ainlO) = limfc_>o ai'i(fc) approaches a constant [60) . Vdd 
connects the incoming partial wave labeled by I with the outgoing partial wave 
labeled by I'. The multi-channel nature of the dipole-dipole potential gives rise 
to novel two-body resonances, which can be tuned by varying the dipole moment 
d [MOl ma [1431 IISII Hm [He]. For electric dipoles, this can be achieved through the 
application of an external electric field. Similarly to the case of purely spherically- 
symmetric two-body interaction potentials, the energy spectrum of two dipoles under 
spherically symmetric harmonic confinement can be determined analytically if the 
dipole length d* is smaller than the harmonic oscillator length aho,/i [115] . This review 
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will not discuss the few-body physics of trapped dipolar systems further. 



3. Hyperspherical coordinate treatment of trapped three-particle system 

This section considers three-particle systems with short-range interactions under 
external spherically symmetric harmonic confinement. For this confinement geometry, 
the center-of-mass degrees of freedom Rem separate off, leaving six relative coordinates. 
The trapped three-body system can be analyzed by a variety of techniques including 
those developed by chemists and nuclear physicists. This section focuses on 
the hyperspherical coordinate approach, which provides an intuitive picture for 
understanding three-body systems with bosonic, fermionic or mixed symmetry. 
Subsection 13.11 introduces the formalism, while Subsecs. 13.21 and 13.31 discuss selected 
aspects of the non-interacting and large scattering length regimes and the Efimov 
regime, respectively. Implications of three-body physics for larger systems are 
discussed in Sees. |4] and [5] 



3.1. General framework 

The hyperspherical coordinate approach transforms the coordinates of the system into 
a single length, the hyperradius i?, and a set of hyperangles, collectively denoted by 
O mi Uni IZB HZl Hi IHi IlSni IISI] ■ For the three-body system in the center of mass 
frame, f2 consists of five hyperangles [152] . The hyperspherical coordinate approach 
has been employed in many contexts such as in the description of autoionizing states 
of the He atom |147] and the Efimov effect |153j . It has also been applied successfully 
to systems consisting of more than three particles [TO 11491 [151] (see, e.g., Subsec. lT^ . 
The hyperradial dynamics can be visualized and interpreted in much the same way as 
that of a diatomic molecule described by a set of coupled BO potential curves. 

To define the hyperspherical coordinates for three-body systems with masses rrii 
and position vectors (i = 1, 2, 3), we introduce Jacobi vectors pi, 

Pi = ri- r2 (20) 

and 

TOlfl -f TOlf2 ^ 

P2 = ■ rg. 21) 

mi + TO2 

The Jacobi vector p^ coincides with i?cin- The hyperradius R, which can be interpreted 
as measuring the overall size of the system, is defined through 

HrR^ ^ fiipl+ P2pI (22) 

where p,R denotes the hyperradial mass (a scaling factor that can be chosen 
conveniently), and where /ii and /i2 denote the reduced masses associated with pi 
and /ojj, respectively, /ii = mim2 / {mi + 1112) and p,2 = (rrii -f 7712)7713/(7711 + m2 + 7713). 
The hyperangle a,0<Q;<7r/2, is defined as 

a — arctan [ 2^^^^ | ^ (23) 

The other four hyperangles are given by the angles {'di,ipi) and (7^2, ¥^2), which 
describe the orientation of the unit vectors pi and p2, respectively. In these 
coordinates, the volume element becomes R^dRd^, where d^l = sin{2a)dadflidil2 
with dOj = sin i) jddj dip j. We note that the definition of the hyperangular coordinates 
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is not unique. The "optimal set of coordinates" depends on the choice of the analytical 
or numerical technique employed. 

In the hyperspherical framework, the relative wave function ip{R, f2) is expanded 
in terms of a set of channel functions ^^{R\il) that depend parametrically on the 
hyperradius R and a set of weight functions F^q (i?) [69l I147j , 



HrA 



5 



vq 



MR^ 

sin(2Q;) 



(24) 



where S denotes the symmetrization operator, which can be written in terms of 
permutation operators Pjk that exchange the j*'' and the fc*'' particle. For three 
identical bosons, e.g., S can be written as 3~^/^(l + P23 + ^3i)- The channel functions 
are eigenfunctions of the fixed- i? hyperangular Hamiltonian i?adia, 
h'sUR}^ 



i?adia$.(i?;f^) 



where 



adii 



Mnt(i?,f1). 



(25) 



(26) 



2^iRR-^ 

For each fixed hyperradius _R, the channel functions form an orthonormal set, i.e., 
/ $*(i?;f1)$i.'(i?;17)dl1 = 5^^' ■ In Eq. denotes the square of the grand 

angular momentum operator, which contains the contributions of the hyperangular 
motion to the kinetic energy operator |145| . For the purpose of this review, the 
functional form of is not needed. The potential energy surface Vint accounts for 
the particle-particle interactions, Vint = X]j<A; ^b('"jfc)- 

The weight functions F^q{R) are solutions to the coupled equations 

Hr - ^W,,{R) - E,} F,q{R) = 



where 



W,,,{R)+2P,,,{R) 



d_ 

dR 



Fu'q{R), 



52 h^[sl{R) -1/4] 1 



2 d2 



fiRUj R 



(27) 



(28) 



2^lR dR^ 2/ii?i?2 

The last term on the right hand side of Eq. ([25)1 represents the external spherically 
symmetric confining potential Vtrap in the relative coordinates, Vtrap = X^i ^JT-iW^ff /2 — 
Mcmw2_R2^y2 with Mem = Yli^i- The coupling elements Pyu'{R) and Wyi,'{R) are 
defined through 



and 



Pu.'{R) 



W,,,{R) = / $.(i?;r!) 



d''< ^,-{R;n) 



dn. 



(29) 



(30) 



No approximations have been made, and Eqs. and (P7|) are equivalent to the 
relative Schrodinger equation. 

Note that T4rap enters into the hyperradial Schrodinger equation, Eq. (P7|) . but 
not into the hyperangular Schrodinger equation, Eq. (|25|) . This implies that the 
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three-body system under spherically symmetric harmonic confinement shares many 
properties with the corresponding free-space system. In particular, the hyperangular 
eigenvalues s^{R) / {2^fjR^) are the same for the free-space and trapped systems. On 
the other hand, the weight functions F^q{R) of the free-space system obey, depending 
on the threshold, either bound state or scattering boundary conditions at large R, 
while all Fuq{R) of the trapped system obey bound state boundary conditions at large 
R. A comparison of the two-body and three-body systems reveals that the three-body 
hyperradius R plays a role analogous to the two-body interparticle distance r, and 
that the hyperangles ^ play a role analogous to the angles that describe the direction 
of r . 

In general, the coupled equations in the hyperradial coordinate R, see Eq. ([27)) . 
have to be solved numerically. However, certain approximations can be made that 
lead to strict upper and lower variational bounds [Ml I154[ 1155] . An upper bound to 
the energies can be obtained by neglecting the coupling matrix elements W^^i and 
Pjyjy' with V ^ v\ i.e., by setting the right hand side of Eq. ([?7)) to zero. The energies 
obtained by solving the radial equation in this approximation, the so-called adiabatic 
approximation, provide an upper bound to the exact eigenenergies. A lower bound to 
the exact ground state energy can be obtained by additionally neglecting the diagonal 
coupling element or so-called adiabatic correction Ww 

If Vint is given by a sum of two-body zero-range potentials, i.e., T^nt = 
Sj<fc ^^''(^ifc)' the eigenvalues fi^s^(i?)/(2^fli?^) of the hyperangular Schrodinger 
equation can be obtained by finding the roots of a transcendental equation |156| for 
any combination of masses mi, any combination of s-wave scattering lengths ao"''^''(0), 
and any symmetry (for earlier work see, e.g., Refs. ^69 l 11531 11571 11581 1159j ). For this 
review, the BBB system (three identical bosons) and the FFX system (two identical 
fermions and a third particle, either boson or fermion) are the most relevant. The FFX 
system with mass ratio k = 1 can be realized by occupying two different hyperfine 
states of ^Li or "^^K. FFX systems with unequal masses of current experimental interest 
include ^Li-^^Sr (k « 0.071), ^Li-S^Rb (k w 0.069), 40K-84Sr {k w 0.48), i^^Yb-^OK 
(k « 4.3), "OR-^Li (k « 6.7), and ^^Sr-^Li (k « 14.5). 

For the BBB system with relative angular momentum i = and s-wave 
scattering length ao(0) [i.e., all cl'^^^ i^i) are identical], the transcendental equation 
reads [1531 UMl USQ] 



r ^ 

Si,(i?)C0S s^{K) — 

^4-y2sin 
ao(0) 



— ;= sm 
V3 



s.(i?)|] . (31) 



Equation pil) is written in terms of yJJT^R as opposed to R since the definition of 
i?^, Eq. p2p. contains the scaling factor /i^j; in other words, [irR? is an invariant. 
Figure [5] shows the quantity Si/(i?) for the BBB system with i = as a function of 
^/i_Ri?/[^/TOao(0)]. Si/(i?) is either purely real or purely imaginary. In particular, the 
lowest Sv(R) value becomes imaginary at Y^/i/ji?/[-ymao(0)] ~ —0.6385. Although Si, 
changes from being purely real to being purely imaginary, the quantity s^, which enters 
into the hyperradial Schrodinger equation [see Eqs. ([27]) and ([281 ]. changes smoothly 
(it simply changes sign). The purely imaginary root has profound implications on the 
system properties. As discussed in more detail in Subsec. [3^ the imaginary root in the 
large |ao(0)| regime gives rise to Efimov physics. Expressions analogous to Eq. (|3ip can 
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Figure 5. (Color online) s^{R) values for the BBB system with L = 0. Dotted 
and solid lines show the purely real and purely imaginary roots s^(R), obtained 
by solving Eq. II31I I. as a function of y//Jfj-R/[v'mao(0)]. 



be obtained for other L values, i.e., there exists a sequence of hyperangular eigenvalues 
s^{R) for each L. Efimov physics is absent for BBB systems with L > 0. 

The FFX system contains two s-wave interactions, one for each of the two FX 
pairs. The implicit eigenequations for Sy{R) read [156j 

s^{R)k^'^ cos [su{R)f\ 



(1 + k)1/2 



(1 



^3/2, 



..(fl)sin-i (tTTt) 



[k(1 + 2k)] 



1/2 



R 



for L 



m ao(0) 
and P51IT77] 

,3/2 



Su[R)\ 



(32) 



K^I^[sl{R)-l]2F, 



3+s„(fl) 
2 ' 



3(1 + k)3/2 
[sliR)-lW/^sm [g.(i?)f] _ 



s.(i?)(l 
IfJ'R R 



■ COS 



-k)1/2 



(33) 



TO ao(0) 

for L = 1, where k denotes the ratio between the mass of the F and X atoms and 
TO denotes the mass of the X atom. Similar expressions can be obtained for L > 1. 
Figure [5] shows the effective potentials s1{R) / {2^piR^) for the FFX system (k = 1) 
with L = Q and 1 as a function of ■y///i^i?/[^/TOao(0)], where ao(0) is positive. For the 
three-body system, a large hyperradius R corresponds to one of two scenarios. Either 
all particles are far apart from each other or one particle is far apart from the other 
two. The former is associated with the atom-atom-atom continuum while the latter 
is associated with the atom-dimer continuum. The hyperradial potential curves that 
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Figure 6. (Color online) Hyperradial potential curves h'^sl{R)/{2fiiiR'^) for the 
FFX system (k = 1) as a function of ^fiiiR/[\/mao{0)] for ao(0) > 0. Dashed and 
solid lines show the potential curves that asymptotically approach the atom-dimer 
threshold for L = and L = 1, respectively. Dotted and dash-dotted lines show 
the potential curves that asymptotically approach the atom- atom- atom threshold 
for L = and L = 1, respectively. As a reference, the thin dash-dash-dotted line 
indicates the atom-dimer threshold. 



approach the atom-dimer threshold are shown by dashed and sohd hnes for L = and 
L — 1, respectively. Figure [S] shows that the L — curve approaches the atom-dimer 
threshold (thin dash-dash-dotted line) faster than the L — I curve. The ordering of 
the hyperspherical potential curves in the small i?/ao(0) regime reflects the values of 
Si, at unitarity [i?/ao(0) = corresponds to the infinite scattering length or unitary 
limit]. As elaborated on further in Subsec. 13.21 (see Fig. [7]), the lowest s^, value for 
L = 1 lies below that for L = at unitarity. The fact that Fig. [5] contains only a 
single atom-dimer threshold is unique to the zero-range potential model employed. 
More realistic two-body potentials support several two-body bound states, and the 
corresponding three-body hyperradial potential curves contain as many atom-dimer 
thresholds as the two-body potential supports bound states. The fact that the lowest 
hyperradial potential curves for L = and L = 1 are purely repulsive for k = 1 
signals the absence of three-body bound states for the FFX system with zero-range 
interactions. 

3.2. Analytical solutions in the non-interacting and large scattering length limits 

This subsection applies the hyperspherical framework to non-interacting three-body 
systems and three-body systems with large scattering length. Throughout this 
subsection, we assume that Vint is given by a sum of two-body potentials Vth{rjk), 

(ik) / \ (ik) 

which are characterized by zero-energy s-wave scattering lengths (0). If all Oq 
are either zero or infinity, the scattering lengths do not define meaningful length scales. 
In the zero-range limit, the hyperradial and hyperangular degrees of freedom decouple, 
i.e., the coupling matrix elements Wi,^' and P^i,' vanish identically [4T | I156| [160) . and 
$y and s„ are independent of R. For finite-range two-body interactions with infinite 
scattering lengths, the decoupling is approximate. In this case, the zero-range solutions 
apply in the regime where R and ao(0) are much greater than tq, where tq denotes 
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the range of Vtb [or, if the Vthirjk) are not the same for the different pairs, the largest 
of the ranges of the Vthirjk)]- 

For now, we assume that the s,^ are known [see, e.g., Eqs. (|5T|) through (|33l) ] and 
independent of R. The hyperradial Schrodinger equation then reads 

' VMR) + i:m^'R' 



F,,{R) = E,gF,,{R), (34) 



where 



Here, q is the hyperradial quantum number [see Eqs. (f27| and ()28p ]. Equations (p4| 
and (I55|) are identical to Eq. (|12l) with Vtb = if the identifications R — r, — fi 
and Sjy — 1/2 = / are made. This subsection assumes that Si, is real and greater 
than zero; the next subsection discusses the case where Sj/ is purely imaginary. The 
two linearly independent solutions of Eq. (j34[) are thus identical to those given in 
Eqs. (|13p and ITU) if the appropriate substitutions are made. The main difference 
between Eq. (IM)) and Eq. (|12p with Vth = is that the quantity 5,^ — 1/2 can take 
non-integer values (s^ — 1/2 > —1/2) while the orbital angular momentum quantum 
number / takes the values 0, 1, • • •. Enforcing the bound state boundary condition at 
large R, the hyperradial wave function is given by Eq. ()16|) with the substitutions 
introduced above. 

To determine the quantization condition for g, we investigate the small R behavior 
of the two linearly independent functions fvq{,Fi) and gvq{R) (here and in the following, 
the first subscript is written as v as opposed to s,y — 1 /2 for notational convenience) . 
Using lim^^o (iFi(a, 6, x^)) = 1 [103] in Eqs. ^ and we find f^q{x) 
and gvq{x) — ?> ^.-Si^+i/z as a; ^ 0, where x = R/a\yo,R with Oho.R = \/f'-/{pR'-^)- fi^q{x) 
is normalizable for all s^, + 1/2 > —1/2 while g^qix) is only normalizable if is less 
than one |161| . The latter follows from the fact that /J \g,yq{x)\'^dx = Iu^Iq + C 
for = 1, where C is an integration constant and e denotes a small positive 
quantity; thus, the integral shows divergent (convergent) behavior in the x ^ 
limit for > 1 (s^ < 1). Correspondingly, the irregular function g^qix) must be 
eliminated for Si, > 1 [ITl 11611 1162) . Enforcing the boundary condition Fi,q{x) — > 
on the exponentially decaying piece of fvq{x), one finds q = 0, 1, • • • for > 1. For 
< Si/ < 1, the irregular solution is normalizable and the hyperradial solution is given 
by a linear combination of the regular and irregular solutions. The ratio between the 
two solutions is determined by the logarithmic derivative or hyperradial boundary 
condition at small R, which can be thought of as being imposed by an effective 
hyperradial three-body potential. If the irregular solution is suppressed naturally, 
which is expected to be the case for "generic interactions" , the properties of the system 
with 1 > Sj/ > are fully determined by ao(0) [16H 1162] . However, as discussed in 
Refs. |163|I1641I165) and Subsec IO] intriguing physics emerges if both the regular and 
irregular solutions contribute. Table [2] summarizes the dependence of the quantization 
condition on the value of s^. 

We now apply the outlined formalism to the BBB and FFX systems |158ill60] . For 
three identical bosons, the wave function has to be fully symmetric under the exchange 
of any pair of bosons. This symmetry constraint reduces the density of states compared 
to that of a XYZ system, which consists of three distinguishable particles that interact 
through the same potential energy surface as the BBB system. Since the hyperradius R 
is unchanged under the exchange of any pair of particles, the symmetrization operator 
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Table 2. Dependence of the solutions to the hyperradial Schrodinger equation 
with /?-indcpondent s^, Eqs. II27I I and II28I I. on the value of s^. 

real, > 1 ij = 0, 1, ■ ■ ■ 

Si/ real, 1 > Sjy > q depends on hyperradial boundary condition 
Su imaginary q depends on hyperradial boundary condition 

(Efimov physics) 



S only affects 0. and, correspondingly, only the hyperangular channel functions. The 
ground state energy of the non-interacting BBB system, which has vanishing relative 
angular momentum and natural parity (i.e., = 0+), corresponds to a situation in 
which all bosons occupy the lowest single particle state. Since each single particle state 
has an energy of 3?iw/2, the relative ground state energy E is E = 3huj (or sq — 2) and 
the center of mass energy is Ecm = 'ihw/2. The allowed Si, values of the excited non- 
interacting states can be found by solving Eq. ([3T|) for ao(0) = (see also asymptotic 
values in Fig. [5]). For L = 0, we find ^ 2, 6, 8, 10, 12, 14, 16, 18, •• • {s„ = 4 is 
forbidden due to symmetry constraints), with degeneracies 1, 1, 1, 1, 1, 2, 1, 2, • • •; the 
Si, values for 18 and higher are doubly-degenerate |166j . This non- interacting L — Q 
example shows that the Su values do not necessarily follow a simple, easy to guess 
pattern. For each Si, (which depends on L), q takes the values g = 0, 1, • • •. 

Next, we consider the BBB system with infinite scattering lengths |158l 1160) . 
The lowest few Sy values for L = are s^ w 1.00624i, 4.46529, 6.81836 and 9.32469 
(see Fig. [5|). To connect the allowed s^, values for vanishing and infinitely large 
scattering lengths, we follow the dotted lines in Fig. [5] from the left to the right. As 
y/jMRR/[^/7TiaQ{Q)\ increases from — oo to —0.6385, the lowest value decreases from 2 
to 0. At y7^^/[v^oio(0)] « —0.6385, the lowest s^, value becomes purely imaginary 
and the dotted line "goes over" to the solid line. The values corresponding to 
excited hyperangular states decrease as ^/ifli?/[-\/TOao(0)] increases from — oo to to 
00. The overall behavior of the Si, values {su real) as a function of ^JJT^R / [^/Triaa{Q)\ is 
similar to that of the eigenenergies of the two-body system as a function of aho,^i/ao(0) 
[see Fig. [3Ka)]. However, the amount by which drops when going from the non- 
interacting to the large scattering length limits depends, unlike in the two-body case, 
on the state under consideration. Furthermore, it must be kept in mind that the 
Si/ values are related to the hyperangular eigenvalues and that coupling between the 
different hyperangular channels modifies the energy spectrum for finite ao(0). 

Next, we consider FFX systems with k = 1 and k > 1. Figure [7] shows the 
Si, coefficients for the FFX system (k = 1) as a function of L for ao(0) = oo 1158) . 
The s^ coefficients are purely real and greater than 1 for all L and correspondingly 
q takes the values 0, 1, • • • for each s,, (see Table [5]). In Fig. [71 dotted lines show 
the approximate semi-classically derived s^, values for oq (0) = oo )158) , = 2i/ -I- 3 
for L — and s^, = 2;/ + L + 1 for L > 1 (;/ = 0, 1, • • •). These expressions, which 
have also been deduced using an "atom-dimer type picture" (1671 1168] , become more 
accurate with increasing v and/or L. The Si, coefficients at unitarity constitute a 
key ingredient for determining the finite temperature equation of state of strongly- 
interacting two-component Fermi gases in the bulk limit via a high temperature virial 
expansion approach (see Subsec. H^ and Refs. [1691 1170) ). 

Figure |8] shows the real part of the lowest Sy coefficients of the FFX system for 
ao(0) = oo as a function of the mass ratio k for L = — 3. The coefficient for 
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Figure 7. (Color online) Symbols show the Su coefficients for the FFX system 
(k = 1) as a function of L for ao(0) = oo. Note that the values are purely 
real for all L. The values have been obtained by solving Eq. (17) of Ref. 11581 
[for L = and 1, see Eqs. I I32II and ||33| |]. Dotted lines show approximate semi- 
classically derived expressions (see text). 




Figure 8. (Color online) Dotted, solid, dashed and dash-dotted lines show the 
lowest so value as a function of k for the FFX system at unitarity for L = 0, 1, 
2 and 3, respectively, and IT = (—1)^. The sq values arc obtained by solving the 
equations given in Sec. Ill of Ref. |156| [for L = and 1, see Eqs. II32I I and Il33l l1 . 



L = 1, which is lower than that for L = 0, 2 and 3, decreases from 1.773 to 1 to 
as K increases from 1 to 8.619 to 13.607. For k > 13.607, the lowest Si, coefficient for 
= 1~ becomes purely imaginary, indicating that Efimov physics plays a role for 
two-component Fermi systems with k > 13.607. Applying the results from Table [5] 
to the FFX system with = 1~ , q takes integer values for k < 8.619, depends on 
the short-range hyperradial boundary condition for 8.619 < k < 13.607, and depends 
on the three-body Efimov parameter for k > 13.607. The fact that the lowest 
coefficient of the FFX system for L = 1 decreases with increasing k can be attributed 
to an effective attraction between the two heavy fermions, which is induced by the 
light particle. This interpretation emerges within a BO treatment |162j . in which the 
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light particle ("electron") moves in an effective potential created by the two heavy 
particles ("nuclei"); thus, the FFX system with large k can be interpreted in much 
the same way as, for example, the molecule. 

3.3. Efimov physics 

As already alluded to, the Efimov effect arises if two or more of the s-wave scattering 
lengths that characterize the two-body subsystems diverge and if Sj^ is purely 
imaginary. Examples include the BBB system with = 0"*" symmetry and the 
FFX system with k > 13.607 and — 1^ symmetry. In the limit that tq — ^ 0, 
the hyperradial potential VeS, Eq. takes the form h'^{—sl — l/4)/(2^fli?^), 

where Si, is purely real, i.e., = lm{s^). The bound state solution of the free- 
space hyperradial Schrodinger, Eq. with w = and Vcs given above, reads 
Fuk{R) — (tiRy^^Kis^^RR), where k — y/2fiii\E\/h and where Kig^ denotes the 
modified Bessel function of the second kind with imaginary index |48| . Here, the 
hyperradial solution is labeled by k and not by k, as done in Subsec. 12.11 since we 
are considering bound states and not scattering states. To be consistent with our 
previously introduced conventions, our R includes an extra factor of \/2 compared 
to Eq. (152) of Ref. [48 . An analysis of the solution shows that the three-body 
system supports an infinite number of three-body bound states. While the bound 
state energies ij'"' depend on the short range boundary condition, the ratio between 
two consecutive bound state energies i?'"' and is independent of the three-body 
parameter and fully determined by Si,, 



For the BBB system in free space, for example, the spacing between the {n+ 1)^^ and 
n}^ bound state is approximately 1/515. This geometric spacing reflects a discrete 
scale invariance and is the key characteristic of the Efimov effect. Other observables 
such as the three-body recombination rate K3 to weakly bound dimers and scattering 
resonances also carry signatures of the discrete scale invariance. For example, K3 can 
be written as /ft|ao(0)|'*/r7i, where the "amplitude" / depends on the two-body s- 
wave scattering length ao(0) and the three-body phase 9b, which is — in turn — related 
to [ml HIS HIS [HI]- The short-range phase 9b can thus be determined by 

fitting the experimentally determined three-body recombination rate, measured over 
a wide range of s-wave scattering lengths [52], to the known functional form of K3. 
Although the discussion so far is, strictly speaking, only valid in the ro — ^ and 
ao(0) — >■ cxD limits, the key features survive as long as ro/|ao(0)| <C 1. See, e.g., 
Refs. [1751 11761 1177] for a discussion of finite-range effects. 

We now investigate how the external spherically symmetric confinement modifies 
the Efimov spectrum |158[ 1160(11781 1179] . In particular, we seek solutions to Eqs. ([M)) 
and p5|) for purely imaginary and i?-independent s^. As before, the hyperradial 
solutions are given by Eq. (fTH]) with I replaced by isi, — 1/2 and x equal to i?/aho.fl. 
The radial solutions can be conveniently written in terms of the Whittaker function 
W [1031 1158] ■ Ft,K.ix) = Nis^{K)x~^/'^WE/2,ig,,/2ix'^), where n is used instead of q to 
label the hyperradial solutions and where Nis^{K) denotes a normalization constant. 
In the small x limit, the radial solutions become 




(36) 



F^i.,{x) N,sjK)^/x{b: 



IX 



b*x-''-) 



(37) 
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Figure 9. (Color online) The solid lines show the relative eigenenergies E, 
calculated using Eq. I)40ll . as a function of di, for the BBB system with = 0"*" 
at unitarity {su 1.00624). For comparison, dashed and dotted lines show the 
eigenenergies for a 10 times smaller and a 10 times larger Sv value, corresponding 
to the FFX system with = 1~ and k just larger than 13.607 and very large, 
respectively. Asterisks mark the {Of,, E) values for which Fig. 1101 shows the radial 
wave functions. 



where [160] 

^-^TiiT^i:^- ^^^^ 

We emphasize that the quantity b depends on the energy of the trapped three-body 
system. Thus, for a given three-body system, b can be varied by changing the angular 
trapping frequency. Defining b = \b\ exp(i6'h), Eq. ((37|) becomes 

F^Kix) Nis^{K)^/xsin{s^ lux + 6b) (39) 

with 

- arg(&). (40) 



For a fixed 9b^ the allowed energy eigenvalues are determined by Eq. pOj) . 

Figure [H] shows the three-body eigenenergies as a function of the three-body phase 
Ob at unitarity for three different Sy values, i.e., = 1.00624, 0.100624, and 10.0624. 
It can be seen that the energy spectrum depends quite strongly on 9b- In particular, 
the energy spacing between consecutive eigenenergies varies notably with Ob for fixed 
Sii. For negative energies, the spectrum looks Efimov-like; by this we mean that 
the spacing between consecutive energy levels follows, very roughly, that of the free- 
space system. For positive energies, in contrast, the spectrum shows similarities with 
the gas-like spectrum of three-body systems for which Efimov physics is absent. We 
emphasize that the three-body phase Ob is not determined by the zero-range model and 
needs to be extracted from experimental data or ab initio calculations. Spectroscopic 
measurements for the trapped three-body system may, e.g., allow one to extract the 
three-body phase 9b |179j . 

The short-range phase 9b determines the small R hyperradial boundary condition. 
The limiting behavior of F^^(x)^ Eq. (P^ . indicates that the radial eigenfunction shows 
oscillatory behavior at small R. Figure [TU] shows the radial solutions -Fi,„ {x) for the 
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Figure 10. (Color online) Solid, dotted, dashed, dash-dotted and dash-dash- 
dotted lines show the radial wave function /2,so/2{^'^)/ for the BBB system 
at unitarity (so = 1.00624) for E = —Ahw, 0, 2hw, Uhw/A and 7tkd/2, respectively. 

BBB system at unitarity {si, — 1.00624) for five different relative energies E/{huj) 
or, alternatively, for five different three-body phases 9^. The corresponding ((?;,, E) 
values are marked by asterisks in Fig. [9] It can be seen that the wave length of the 
oscillations decreases with decreasing x and that the wave function extends to larger 
hyperradii with increasing energy. 

A more detailed discussion of Efimov physics can be found in Refs. |48j|49]. The 
key point for the discussion of larger systems is that the three-body energy spectrum 
in the Efimov regime depends, in addition to the scattering length, on a three-body 
parameter (parameterized through 9b above). Although the behavior at unitarity is 
particularly transparent since the channel coupling vanishes, the three-body phase also 
plays a role for systems with finite scattering lengths. As shown in Subsec. 15.31 the 
energy of weakly-interacting trapped Bose gases also depends — although in a rather 
different manner — on the three-body parameter. 

4. Fermi gases under external confinement 

This section discusses the behavior of s-wave interacting Fermi gases under external 
confinement. Subsection 14.11 provides a qualitative and, where possible, quantitative 
description of the crossover from the weakly- attractive atomic Fermi gas to the 
weakly-repulsive molecular Bose gas regime of equal-mass systems and introduces 
a perturbative description of the weakly-interacting regimes. Subsection 14.21 focuses 
on the unitary regime, where the s-wave scattering length diverges. Subsection 14.31 
extends the discussion of equal-mass Fermi gases to unequal-mass Fermi gases. Lastly, 
Subsec. Hm shows that the microscopic few-body energy spectra can be used to predict 
the thermodynamic properties of inhomogeneous and homogeneous equal-mass two- 
component Fermi gases via a high temperature virial expansion approach, which — 
somewhat surprisingly — remains valid down to about half the Fermi temperature. 
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Microscopic description of BCS-BEC crossover of trapped two-component Fermi 
gas with equal masses 

This subsection considers two-component Fermi gases with Ni atoms of the first 
species and N2 atoms of the second species {N = Ni -\- N2). We assume that like 
fermions are non-interacting and that unhke fermions interact through a sphericaUy 
symmetric short-range two-body potential Vth{rjk), where rjk — \rj — rk\ [rj denotes 
the position vector of the j"* fermion (j = 1, • • • , N) measured with respect to the trap 
center]. Restricting ourselves to spherically symmetric harmonic confining potentials 
with angular trapping frequency uj, the model Hamiltonian H reads 



where m denotes the atom mass. The Hamiltonian given in Eq. (|^T|) describes equal- 
mass two-component Fermi gases such as ®Li atoms in two different hyperfine states 
or '^^K atoms in two different hyperfine states. The Hamiltonian is single- channel 
in nature, which — following the discussion of Subsec. 12.11 — implies that our treatment 
applies to Fermi gases near broad but not near narrow Fano-Feshbach resonances. This 
subsection discusses the properties of systems governed by H, Eq. (PT|) . as functions 
of the number of particles N (throughout, we use Ni > N2) and the s-wave scattering 
length ao(0) of the two-body potential Vtb- 

Equal-mass two-component Fermi gases are fully universal, i.e., the system 
behavior is fully determined by ao(0), provided the range rg of the two-body potential 
Vtb is much smaller than the s-wave scattering length ao(0), the average interparticle 
spacing (r^j) and the harmonic oscillator length aho, where Oho = \/h/ {muj). The 
evidence for universality comes from studies that show that the system behavior 
remains essentially unchanged when the shape of the two-body interaction potential is 
varied in the regime rg <C Oho and tq ^ {fij)- This characteristic is closely linked to the 
fact that the (2,1) system does not support a weakly-bound three-body state |162[ll80j 
(see Subsecs. 13.11 and 13.21) . Furthermore, weakly-bound four-body or higher-body 
bound states are absent in the zero-range limit [37l |38l I181i I182i 1183] , and dilute two- 
component Fermi gases are stable even for infinitely large s-wave scattering length 
The Pauli exclusion principle can be thought of as producing an effective repulsive 
force that stabilizes the system with attractive interactions against collapse. This is 
reminiscent of the stabilization of white dwarfs against gravitational collapse by the 
electron degeneracy pressure |184j . 

We start our analysis of the Hamiltonian H, Eq. (|4T|). by considering the 
{Ni,N2) — (2, 1) system. The energy spectrum of the (2, 1) system can be obtained 
using the hyperspherical framework outlined in Subsec. 13.11 Here, we instead 
employ an approach based on the Lippmann-Schwinger equation, which allows for 
the determination of the energy spectrum of the equal-mass three-fermion system 
with zero-range interactions with comparatively little computational effort |185) . The 
success of this technique rests in the fact that two-body correlations are build into 
the three-body wave function from the outset. Figure [TT] shows the energy spectra for 
the two lowest relative angular momenta, i.e., for L — and 1, of the three-fermion 
system as a function of aho/ao(0). In determining the eigenenergies, we assumed that 
the center-of-mass degrees of freedom are in the ground state, i.e., that i?cm = 3hw/2. 
The key characteristics of the three-fermion spectra can be summarized as follows: 

(i) In the ao(0) 0^ limit, the three-fermion system behaves like a non- 




(41) 
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Figure 11. (Color online) Energy spectra for the FFX system with ft = 1 for (a) 
= O'^ and (b) = 1~ as a function of the inverse scattering length [ao(0)]~^ 
{Etot includes the center-of-mass energy of 3/kj/2). The BCS-BEC crossover curve 
is shown by dashed lines. The energy spectra are calculated numerically following 
the formulation introduced in Ref. |185| . 

interacting atomic Fermi gas. The energy of the lowest L = and 1 states is 
E^at = 13fia;/2 and llSa;/2, respectively. The fact that the lowest L = \ state has a 
lower energy than the lowest L = state can be understood readily by considering 
the non-interacting limit, where the wave function separates into a component that 
depends on ri2 and a component that depends on the second Jacobi vector (when the 
interactions are turned on, the wave functions no longer separate and the argument 
needs to be modified accordingly). The three- fermion wave function is antisymmetric 
under the exchange of the two spin-up atoms, implying that the wave function 
component along the ri2 vector carries one unit of angular momentum. For the lowest 
state with L = 1, the wave function component along the second Jacobi vector carries 
no angular momentum (implying that the wave function component along the second 
Jacobi vector contributes an energy of 3hu!/2). For the lowest state with L = 0, in 
contrast, the wave function component along the second Jacobi vector carries one unit 
of angular momentum, which couples to the angular momentum associated with the 
ri2 coordinate such that L = (implying that the wave function component along the 
second Jacobi vector contributes an energy of 5huj/2). 
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(ii) In the ao(0) — !■ 0+ limit, the three-fermion spectra consist of two different 
energy famihes. The first energy family "dives down" as ao(0) approaches 0"'" 
(this is the "molecule+atom family") while the second energy family, the so-called 
"atom+atom+atom family" , approaches the energy spectrum of the non-interacting 
atomic Fermi gas [i.e., the atom-|-atom-|-atom part of the spectrum is the same as 
that in the ao(0) — > 0~ limit (see above)]. The lowest energy level of the (2, 1) system 
belongs to the molecule-|-atom family, for which the system consists of a tightly-bound 
s-wave molecule and a spare atom that carries the angular momentum L of the three- 
body system. The energies of states belonging to the molecule-|-atom family approach 
Etot — -Abound + (2neff + L-\- 3)hu!, whcrc Ucs = 0, 1, • • • and where inbound denotes the 
energy of the tightly-bound s-wave molecule [see Eq. ([7])], as ao(0) — > 0+. Thus, the 
lowest L = state has an energy that is approximately hoj below that of the lowest 
L ~ 1 state. The energy shift due to the effective interaction between the molecule 
and the atom is discussed below [see Eq. dS])]. 

(iii) As discussed in (i) and (ii), the lowest energy of the (2, 1) system has angular 
momentum L = 1 in the ao(0) — > 0~ limit and L = in the ao(0) — )• 0+ limit. 
The lowest L = 1 and L = states cross at aho/ao(0) « 1 [HSl HHSl [M]. The 
energy branch shown by dashed lines in Figs. [TlT a) and [TlT b) corresponds to the 
trap analog of the BCS-BEC crossover curve. While the BCS-BEC crossover curve of 
the homogeneous system depends on two dimensionless parameters, i.e., kpao{0) (kp 
denotes the Fermi vector) and the population imbalance p, p = {Ni — N2)/N, that 
of the harmonically trapped system depends on three dimensionless parameters, i.e., 
oo(0)/aho, N and iVi — N2- The system properties of the (4, 2) and (6, 3) systems, for 
example, differ somewhat, despite the fact that both are characterized hy p = 1/3. 

(iv) The (2, 1) energy spectra show sequences of avoided crossings in the strongly- 
interacting regime, i.e., for aho/|ao(0)| ^ 1. These avoided crossings have been 
analyzed within a Landau-Zener framework [188j and play a crucial role in time- 
dependent studies. One can, for example, prepare an initial state at equilibrium and 
analyze the system response to slow (adiabatic) and fast (non-adiabatic) magnetic 
field ramps, modeled using a time-dependent s-wave scattering length. At the end of 
the field ramp, the occupation of the energy families, and branches within the energy 
families, can be monitored. Calculations along these lines are of direct relevance 
to ongoing experiments. On the one hand, these calculations provide insights into 
the characteristic time scales that may be extrapolated to larger systems. On the 
other hand, these calculations describe the dynamics of an isolated optical lattice site. 
Knowledge of the single optical lattice site dynamics constitutes a crucial building 
block for understanding the dynamics of systems in which lattice sites are coupled. 

The features pointed out in (i)-(iv) for the (2, 1) equal-mass system apply, with 
appropriate modifications, also to larger trapped two-component Fermi gases. In 
particular, for small |ao(0)|, ao(0) < 0, the gound state of the trapped (iVi, 7V2) system 
behaves like a weakly-attractive atomic Fermi gas. For small |ao(0)|, ao(0) > 0, the 
ground state of the trapped system behaves like a weakly-repulsive molecular Bose 
gas consisting of N2 diatomic s-wave molecules and A^i — N2 impurity fermions (if 
A^i ^ N2, it is more appropriate to refer to the molecules as impurities immersed in 
a Fermi sea). In both regimes [negative and positive ao(0)], perturbative treatments 
reveal the leading order behavior of the energy levels. 

For small ao(0), the unperturbed eigenstates of the iV-atom family are given 
by the non-interacting gas-like atomic states and the sum of atom-atom potentials 
between unlike fermions is treated as the perturbation. The non-interacting wave 
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Figure 12. (Color online) Coefficients c{Ni,N2) that determine the energy shift 
of the energetically lowest lying state of the weakly- interacting Af-atom family 
[see Eq. Il43l l] as a function of N. Circles and squares correspond to A^i — N2 = 
and 1, respectively; the solid line serves as a guide to the eye. The values of the 
coefficients are taken from Table 1 of Ref . I187| . 



functions can be written as a product of two Slater determinants (one for each 
component), which are constructed from the single particle harmonic oscillator 
orbitals |187l I189| . If the atom-atom interaction potential is approximated by the 
Fermi pseudopotential Vp with s-wave scattering length ao(0) [78j , 

Vpirjk) = S{rjk), (42) 

where fj, denotes the reduced mass of the two colliding bodies, the leading order energy 
shifts can be calculated semi-analytically within first order degenerate perturbation 
theory. Note that Vp is identical to Vq^, Eq. (j4]), without the regularization operator 
{d/drjk)rjk. Diagonalization of the perturbation matrix results in the approximate 
expression for Etot{Ni,N2) [I871IT89] . 

Etot{m,N2) - <U^i, ^2) + ciN,,N,)^hLo, (43) 

Oho 

where E^J^.{Nl, N2) denotes the total energy of the unperturbed (non-interacting) 
atomic state under consideration and c{Ni,N2) a coefficient. Figure [T^ shows the 
coefficients c{Ni,N2) corresponding to the lowest energy state for systems with iV < 20 
and iVi — A^2 = or 1. Interestingly, the c{Ni,N2) coefficients show distinct even-odd 
oscillations, with c being larger for even N than for odd N. If ao(0) is negative, as in 
the weakly-interacting BCS regime, then a larger c coefficient corresponds to a lower 
energy. Thus, the perturbative energy expression predicts the existence of a finite 
energy gap AE{N) in the weakly-interacting regime [see Sec. 14.21 Eq. for the 

definition of AE{N)]. 

For small and positive ao(0), the perturbative description of the lowest energy 
family is based on the assumption that N2 s-wave molecules, consisting of one spin-up 
and one spin-down atom, form and that the remaining iVi — N2 spin-up fermions are 
unpaired. Effectively, the Fermi gas consists of two constituents (fermions of mass m 
and molecules of mass 2m) and the system properties are governed by the molecule- 
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Table 3. Summary of selected studies that determined the three-dimensional 
atom-dimer and dimer-dimer scattering lengths a^d and a^^, respectively. 



aad/ao(0) 


method 


K = 1: aad/ao(0) 


= 1.2 


free-space scattering 180j 


any k 




free-space scattering 162] 


K = 1: aad/ao(0) 


= 1.1790662349 


free-space scattering 42 


K = 1: aad/ao(0) 


= 1.18(1) 


trap spectrum 187 


add/ao(0) 


method 


K = 1: add/ao(0) 


= 0.6^'') 


free space scattering (four-body) [182) 


K = 1: add/ao(0) 


= 0.60(') 


free space scattering [190J 


K = 1: add/ao(0) 


= 0.60('') 


diagramatically 119111192] 


1 < K < 13.607 




free-space scattering |99| 


K = 1: add/ao(0) 


= 0.608(2) 


trap spectrum [TsT', 189] 


1 < K < 20 




trap spectrum 187 189 


K = 1: add/ao(0) 


= 0.62(1) 


homogeneous system 38' 


K = 1: add/ao(0) 


= 0.605(5) 


free-space scattering 193 



(<»)The accuracy is reported to be 2 %. (''^No errorbar is reported. 



molecule (or dimer-dimer) scattering length add and the atom-molecule (or atom- 
dimer) scattering length Oad- Following this logic, the unperturbed wave function is 
written as a product of a permanent, constructed from N2 harmonic oscillator orbitals 
with width yjh/{2muj) that depend on the center of mass vectors of the molecules, 
and a Slater determinant, constructed from N1 — N2 single particle harmonic oscillator 
orbitals with width ^Jh/ [muj). The interaction between an atom and a molecule 
(between two molecules) is parameterized by Eq. (|42|) with ao(0) replaced by aad 
(add) and /i replaced by /Xad = 2m/3 {^dd = w). For the lowest energy family, the 
perturbative treatment gives [1871 1189] 

^tot(A^l, A^2) « A^2i^dimcr + <U^1 ^ A^2, 0) + 
N2(N2-l) [2 Odd , . flad j_ 

\ hi^ + N2(Ni - N2)Ci,d nw, (44) 

^ V TT aho,^dd «ho,Mad 

where aho,fii = \/hJ{piUJ) with i = ad or dd. In Eq. pi)) . i^dimor denotes the ground 
state energy of the trapped dimer interacting through Vq^ with scattering length ao (0) 
[i.e., the lowest / — eigenenergy of Eq. (IT7)) . with the center of mass energy of 3?iaj/2 
added]. The coefficient Cad depends on how many unpaired fermions there are and 
which single-particle states the unpaired fermions occupy. For the ground state of 
systems with A^i — N2 = and 1, we have Cad = and Cad = \/2/7r, respectively. 

For the (2, 2) system, for example, Eq. describes the energies of the 

molecule-t- molecule family [ao(0) small and positive] while Eq. (i43l) describes those of 
the four-atom family [|ao(0)| small and either positive or negative]. The (2, 2) system 
additionally contains a molecule-f atom+atom family, whose limiting ao(0) — >■ 0"*" 
behavior can be obtained by assuming that only one dimer forms. 

Table [3| summarizes the values obtained for aad and add in the literature by a 
variety of means. Here, we discuss the determination of the dimer-dimer scattering 
length in more detail. In the first Born approximation, the dimer-dimer scattering 
length add equals 2ao(0). Microscopic free-space scattering calculations confirm that 
add is directly proportional to the atom-atom scattering length and that no other 
parameter enters but find that the proportionality factor is ~ 0.6 as opposed to 
2 [182j . The microscopically derived proportionality factor has, for example, been 
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confirmed using a diagrammatic approach |191[ I192[ and by analyzing the equation 
of state of the homogeneous two-component Fermi gas, calculated using Monte Carlo 
techniques |38| , within a perturbative framework analogous to that given in Eq. (j44p . 
Here, we show that the dimer-dimer scattering length can alternatively be extracted 
from the energies of states corresponding to the dimer-dimer family of the trapped 
(2,2) system [l87l[T89] . 

We start our discussion with the normalized energy crossover curve Ani N2 |187l 



EtotiNi,N2 



3{Ni-N2)nw 
2 



^"^^"^ ^ E-^[N,,N,)^lNhu. ' 

which depends on A^i, N2 and ao(0). For the BCS-BEC branch, kj^^^j^^ changes from 
1 to as l/ao(0) changes from —00 to 00. In addition to providing a convenient 
means to plot the energies throughout the crossover, A^-^^jf^ significantly reduces the 
dependence on the underlying two-body potential [ 1871 11891 1194] . 

The solid line in Fig. [T51 shows A2.2 as a function of 1/0^(0). For as{0) — > 0^, 
Etot{Ni, N2) equals E^J^^{Ni, N2) and E'dimcr equals 3huj (we have two fermions with 
energy 3haj/2 each), yielding A2.2 = 1- For 0^(0) — > 0"'', the total energy can be 
thought of as being due to N2 dimers (with relative energy 3fia;/2 and center of mass 
energy 3fiw/2), yielding A2^2 = 0. Inserting Eq. (H^ into Eq. (P5|) . we find that 
A2,2 « (27r)-i/2add/ah o,Aidd i'^ the ao(0) ~^ 0"*" regime. Thus, add can be extracted 
from the full four-particle energy spectrum by fitting the crossover curve to the limiting 
expression. One arrives at the same result by expanding the eigenenergy of the lowest 
energy branch of Eq. pT|) with I — and /i and aho.^i replaced by /idd and aho.^idd > 
respectively, around ao(0) — 0+ (in this approach, the internal dimer energy and 
the center-of-mass energy need to be added). The analysis can be refined by using 
the full two-body energy spectrum as opposed to the Taylor-expanded expression 
and by allowing add to be energy-dependent via Eq. (|3]). The refined analysis gives 
add = 0.608(2)ao(0) |187| 1189] in good agreement with the microscopic free-space 
results (see Table [3]) and also provides an estimate of the dimer-dimer effective range 
'"eff.dd; ''efi,dd = 0.13(2)ao(0) [187, 189] . The positive value of Tcff ,dd is consistent with 
the conclusion by Petrov et al. |182| that the two dimers interact through an effective 
potential with soft core repulsion. Dotted and dashed lines show A2,2 obtained using 
the perturbative energy expressions, i.e., Eqs. (|43| and p4|) . It can be seen that the 
perturbative treatment provides a fairly accurate description for aiio/ao(0) < —4 and 
aho/ao(0) > 1.5, respectively. 

The energy spectra and crossover curve shown in Figs. [TT] and [T3] are for 
the (2, 1) and (2, 2) systems. Crossover studies for systems with > 4 can be 
found in Refs. j l891 11941 1196] . The energies of the (2,2) system were obtained 
using the stochastic variational approach |197j . a basis set expansion approach that 
results in stable and highly accurate results for trapped systems with up to A^ = 6 
fermions throughout the entire crossover regime [IIZl [HSl [Wl [1881 EHSl [IMl [T98] . 
This numerical approach cannot, at least at present, be applied to significantly 
larger systems since the computational effort increases notably with increasing 
N . Alternative numerical microscopic approaches include an effective field theory 
approach [1861 11991 1200] . an effective interaction approach |201j . and Monte Carlo 
approaches. The applicability of the effective field theory approach is presently also 
limited to relatively small particle numbers (A'^ < 4) while Monte Carlo methods have 
been applied to larger systems. The Monte Carlo techniques can, for the purpose of 
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Figure 13. (Color online) Normalized energy crossover curve A2,2 for k = 1 
as a function of aho/(io(0). The solid line shows A2,2 derived from the L = 
extrapolated zero-range ground state energies reported in the supplemental 
material of Ref. 1167). The dotted and dashed lines are obtained using the 
perturbative expressions [Eqs. II43I I and 1144 1 11 in Eq. I I45II . 



this review, be divided into two classes. The first class, which includes the diffusion 
Monte Carlo (DMC, sometimes also referred to as Green's function Monte Carlo) 
and variational Monte Carlo (VMC) approaches, is applicable to essentially any 
interaction strength pMl [Wl [TMl [TMl . These techniques employ the fixed node 
approximation and result in variational upper bounds to the exact eigenenergies |195| . 
Note that the accuracy may diminish for certain regions of interaction strengths. 
The second class of Monte Carlo approaches such as the lattice MC approach from 
Refs. [2021 1203] have thus far only been applied to the unitary regime. 

4-2. Large scattering length regime: Unitarity 

This subsection considers two-component equal-mass Fermi gases with infinitely large 
s-wave scattering length under external spherically-symmetric harmonic confinement. 
As discussed in Subsec. 13.21 an infinitely large s-wave scattering length does not 
define a meaningful length scale for the system. Thus, two-component Fermi gases 
with infinite ao(0) and zero-range interactions are characterized by the same number 
of length scales as the non-interacting system. This implies, as in the three-body 
case, that the hyperradial and hyperangular degrees of freedom separate in the limit 
of zero-range interactions with infinitely large ao(0) [4TJ I204j . Correspondingly, the 
hyperradial Schrodinger equation is given by Eq. p4p with R and fiR now denoting the 
hyperradius and the hyperradial mass of the TV-body system (i.e., fiRR^ — f^iPf^ 
where the pi denote the Jacobi vectors — excluding the center of mass vector — and the 
Hi denote the associated Jacobi masses). Thus, the hyperradial solutions discussed in 
Subsec. [3.21 applv to two-component Fermi gases with arbitary {Ni,N2). For equal 
mass systems with zero-range interactions, the Si, are all greater than 1. This implies 
that the energies i?"ot'* unitarity can be written as [ITJ 1204] 

E^^;' = {2q + s, + l)hiu + E,^, (46) 
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where the hyperradial quantum number q takes the values 0,1,--- and where the 
Si,, which depend on Ni and N2, are determined from the solution of the 3A^ — 4 
dimensional hyperangular differential equation. The separability of the hyperradial 
and hyperangular degrees of freedom implies that the dynamics of the TV-body system 
is governed by a set of uncoupled effective one-dimensional hyperradial potential 
curves. Thus, the strongly-interacting A'^-body problem reduces to an effective one- 
dimensional problem. 

Equation (|46| predicts the existence of energy ladders with spacing 2hu) [HI 1204) . 
Correspondingly, the energy spectrum contains breathing mode frequencies at integer 
multiples of 2fvjj, in addition to frequencies that are associated with transitions 
between different hyperradial potential curves. The latter class of frequencies depends 
on — s^i , since frequencies in this class arise due to transitions between two different 
hyperradial potential curves. These predictions are a direct consequence of the scale- 
invariance [101 El 1204] and can be verified experimentally. 

The characterization of strongly-interacting Fermi gases is non-trivial. Various 
approaches, including those discussed at the end of the previous subsection, have 
been applied. Unlike in the three-body case, for which analytical solutions to 
the hyperangular Schrodinger equation are known (see Subsec. 13. ip . no analytical 
solutions are known for larger systems and one generally has to resort to solving the 
hyperangular equation or the full Schrodinger equation numerically. For systems with 
L = and up to iV = 4 particles, the hyperangular Schrodinger equation has been 
solved by adopting a hyperspherical correlated Gaussian approach [205j . a variant of 
the stochastic variational approach discussed earlier. For larger systems, in contrast, 
obtaining the solution to the hyperangular Schrodinger equation (which would result 
in the s^) is more challenging than obtaining that to the full Schrodinger equation 
(which results in i?"^"'*). Using Eq. p6)) . which has been proven to hold if the range 
ro of the underlying two-body potential is sufficiently small [1831 11981 1206] , i?"oT* 
can be converted into each other. 

Figure [Hf a) shows the total energy E^^^^ as a function of N for two-component 
Fermi gases with TVi — iV2 = or 1 at unitarity, calculated using three different 
approaches (i.e., a fixed node DMC approach [183] . a lattice MC approach [203] 
and a density functional theory (DFT) approach [207| ). These sets of energies have 
been chosen to illustrate our current understanding (see Refs. [1941 1208] I209j for 
other calculations). The MC and DFT approaches are based on entirely different 
formulations and the relatively good agreement displayed in Fig. 114( a) is encouraging. 
To interpret Fig. fWa'). we relate the energies i?"o"'' at unitarity to the energies 
of the trapped non-interacting system. Since the system at unitarity is characterized 
by the same number of length scales as the non-interacting system, the energies have 
to be directly proportional to each other [5U] . 



the reason for writing the proportionality factor as a square root will become clear 
below. The total energy E'^}. of the trapped non-interacting system shows shell 
closures for Ni = 1, 4, 10, 20, • • •, i.e., the energy increases by 3?iw/2 when the first 
spin- up fermion and spin-down fermion are added, by 5fta;/2 when the second through 
fourth spin-up and spin-down fermions are added, and so on. The energies at unitarity 
show — on the scale chosen in Fig. lMT a) — no signature of shell closures but instead odd- 
even oscillations. Thus, it appears reasonable to determine the proportionality factor 
by "smoothing" the non-interacting energies, i.e., by approximating Ef^^ by the 




(47) 
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Figure 14. (Color online) Energetics of trapped two-component Fermi gas with 
AC = 1 at unitarity as a function of N for A'^i — N2 = or 1. Panel (a) shows the 
total energy i?"^'*. Panels (b) and (c) show the energy difference E^JJ'* — i?flt 
using Eflt = v0.467-E(.jjj' (as a guide to the eye, dotted lines connect the data 
points for even and odd A'^). Circles, squares, and diamonds show results based 
on the fixed node DMC energies from Ref. [183| . the lattice MC energies from 
Ref . |203l , and the DFT energies from the supplemental material of Ref. I207| . 



semi-classical extended Thomas Fermi energies E'^J.'^'^^ [210] 

pNI,ETF 



E 



hio, (48) 

where cetf — 1 (in Subsec. 14.31 we treat cetf as a fitting parameter). We refer to 
the energy obtained by fitting the even N energies to Eq. dH]) with E^J^ replaced by 
E'tot^'^'" as Efit. For the energies for iV = 2 - 30 of Ref. [US], for example, ^tr is 



In the large iV limit, the term proportional to cetf can be 
reduces to the local density approximation (LDA) expression 



found to equal 0.465. 
neglected and Eq. 
(see below) [39j. 

Symbols in Figs. [Ml b) and (c) show the quantity 



unit 
tot 



Efjt. The energy 
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difference — -Efit reveals clear odd-even oscillations for all three data sets. For 

the fixed node DMC data [circles in Fig. [HTb)]. the energy difference i^™* ~ ^f^t 
is approximately equal to zero for even TV and of the order of hto (with a slight 
overall increase with increasing N) for odd N . For the lattice MC data set [squares in 
Fig. [TlT b)] , in contrast, i?"o"' ~ ^fit shows oscillations that follow the shell structure of 
the non-interacting system for both even and odd N (recall, the non-interacting system 
with A^i = N2 possesses shell closures at = 2, 8, 20, • • •), with a roughly constant 
offset between the even and odd N data. Thirdly, the DFT data [see Fig. HDJc)] show 
structure that appears to be independent of the shell structure of the non-interacting 
system. It is currently not fully clear which of these behaviors can be attributed to 
numerical artifacts and which reflect genuine physics. We speculate that the data from 
Ref. [203j overestimate the shell structure but note that further studies are needed 
to fully understand the energetics of trapped two-component systems at unitarity. 
We note that modified versions of Eqs. (|T7)) and (pS)) . which are based on a different 
interpretation of the underlying physics, lead to somewhat different ^tr values [211j . 

To further interpret the results of the trapped system at unitarity, we relate 
the energies of the trapped and homogeneous systems via the LDA. The energy 
-^hom P^^ particle of the homogeneous system at unitarity is directly proportional 
to the energy Eyg of the non- interacting homogeneous gas, iJJj™* = ^hom-EpG- The 
proportionality constant ^hom, also referred to as the Bertsch parameter |212| . has 
been determined theoretically and experimentally. The most reliable estimate comes 
from very recent fixed node DMC calculations, which determine an upper bound for 
Chom, 610m = 0.383(1) [213] . The value of ^hom is of fundamental importance and 
enters into effective theories that use the equation of state of the homogeneous system 
as input. Application of the LDA (see, e.g., Ref. [214] ). which treats the trapped Fermi 
gas as being locally (i.e., at each distance from the trap center) characterized by a 
constant density, predicts that the energy of the trapped system at unitary is given 
by i?"o"* = \/6iom£^tot [22] ■ Thus, ^tr (see above) provides an alternative estimate of 
^hom- The fact that the ^tr obtained by analyzing the energies of the trapped system 
(^tr = 0.465) is larger than the value of ^hom = 0.383(1) [213j obtained from recent 
MC simulations for the homogeneous system may be attributed to the fact that the 
calculations for the trapped system are restricted to relatively small and/or that 
the energies of the trapped system may be slightly too high. The latter may be a 
consequence of the fixed node approximation employed in Ref. |183j . 

Another quantity of interest is the excitation gap /S.E{N), which characterizes 
the odd-even oscillations of two-component Fermi gases. The excitation gap combines 
the energies of the imbalanced system with Ni — N2 = \ (i.e., N = 2Ni — 1) with the 
energies of the next smaller and next larger balanced systems [39], 




Figure [15] shows AE{N) as a function of N for two-component Fermi gases at 
unitarity for the same data sets as shown in Figs. [T4][a), [T4lb) and [TlT c). While 
the energy difference -E™' — Efn shows a distinct dependence on the numerical 
approach employed, the excitation gap determined by the three different methods 
agrees fairly well. The excitation gap AE{N) increases from about O.lHuj to about 
l.bhuj as N increases from 3 to 29. In the homogeneous system, the excitation gap 
equals half the energy it takes to brake a pair |39) . In the trapped system, the 




AEiN)=Etot{Ni,Ni-l)- 
EtotjNi - 1, iVi - 1) + EtotjNuNi) 
2 



(49) 
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Figure 15. (Color online) Excitation gap AE{N) of trapped two-component 
Fermi gas with k = 1 at unitarity as a function of A'^ for A^i — N2 = or 1. Circles, 
squares, and diamonds show results based on the fixed node DMC energies from 
Ref. |183) . the lattice MC energies from Ref. |203| . and the DFT energies from 
the supplemental material of Ref. 12071 . For clarity, errorbars are only shown for 
the fixed node DMC and lattice MC data sets. 



interpretation is not quite as straightforward because of the inhomogeneity induced by 
the confining potential. In particular, the fixed node DMC calculations suggest that 
the extra particle is located near the edge of the cloud as opposed to the center of the 
cloud [183[ 1187] . This finding supports the conclusion put forward by Son |215| that 
AE{N) scales as N^^^ and not as N^^^ as would be expected based on the LDA. The 
microscopic data available to date do not, however, extend to sufficiently large N to 
unambiguously distinguish between the N^^^ and N^^^ behaviors. We note that the 
definition of the excitation gap AE{N), Eq. P^ . applies to all interaction strengths, 
and not just the unitary regime. The excitation gap's behavior in the weakly-attractive 
regime has been discussed briefly in the context of Fig. [121 

Lastly, we use the energies i?™' to arrive at an interpretation of the unitary 
Fermi gas within the hyperspherical framework. The odd-even oscillations of the 
energies imply odd-even oscillations of the Si, coefficients that describe the ground 
state and thus an odd-even staggering of the hyperspherical potential curves Vcft{R) 
that govern the dynamics at unitarity [41] 1187] . Application of the LDA predicts that 
the Su coefficient of the ground state for even N approaches V^iom-Bj^}./ (huj) in the 
large N Hmit [183]. 

The microscopic energies as well as structural properties have been used to 
assess the accuracy of a number of effective approaches, such as those based on 
fractional exclusion statistics [216], DFT [2071 HOSl HH HH [220] and effective 
field theory. Lastly, we note that effort has also been devoted to understanding 
population-imbalanced Fermi gases with A'^i ^ N2 at unitarity within a microscopic 
framework \221\ . In the limit of "extreme" population imbalance, the trapped analog 
of the fermionic polaron problem is realized |2221 1223) . 

4-3. Mass-imbalanced Fermi gas 

This subsection considers s-wave interacting two-component Fermi gases with unequal 
masses. While equal-mass Fermi gases are, by now, relatively well understood, our 
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Figure 16. (Color online) Dimer-dimer scattering length a^^ and effective range 
''ofl.dd as a function of the mass ratio k. Solid lines show ajd obtained by solving 
the free-space Schrodinger equation with zero-range interactions (the results are 
taken from Fig. 3 of Ref. 1991 ) while dotted lines show a^d obtained by applying 
the BO approximation [see Eq. (14) of Ref. |229) ] . Circles and squares show ajd 
and r^fj (jd extracted from the energy spectrum of the trapped (2, 2) system (the 
values of and refl dd are taken from Table II of Ref. I187| ) . 



understanding of unequal-mass systems is much less complete. One of the reasons 
is that experiments on unequal-mass systems are just now becoming available |2241 
12251 12261 12271 1228) . Dual species experiments can be hampered by large three-body 
loss rates, which reduce the system's lifetime and thus make it challenging to cool 
dual-species systems to degeneracy. Trapped dual-species systems with mass ratio k 
(k = mi/m2, where denotes the atom mass of the i*'' component) are, in general, 
characterized by two different angular trapping frequencies, i.e., wi ^ loi. If iO\ and 
iOi differ, the center of mass motion no longer separates. The behavior of the (2, 1) 
system changes with increasing k (see Sec. [3]): Three-body resonances may occur for 
8.619 < K < 13.607 gll HH] and Efimov physics plays a role for k > 13.607 051112]. 
The treatment of larger systems with k > 8.619 thus requires theoretical tools that 
are capable of accounting for these three-body phenomena and, at the same time, of 
capturing the new physics introduced by the fourth, fifth, etc. particle. 

We first consider unequal-mass Fermi gases with small, positive ao(0). Assuming 
that ao(0) fully characterizes the system dynamics, the perturbative expression 
Eq. (|44)) may be generalized to unequal- mass systems |121[ I189j . Thus, a natural 
question is how the atom-dimer and dimer-dimer scattering lengths Cad and add vary 
with K. For the (2,1) system, the atom-dimer s-wave scattering length is associated 
with the = 0+ wave function; for this channel, Efimov physics is absent. The atom- 
dimer scattering length Oad equals 1.18ao(0) for k = 1 and increases with increasing 
K [162] . For the (2,2) system, the s-wave dimer-dimer scattering length was first 
determined via a free-space calculation that employs a zero-range interaction model 
and assumes the absence of three-body bound states for k < 13.607 [99]. add changes 
monotonically from 0.6ao(0) to 1.14ao(0) as k increases from 1 to 13.607 [99^. The 
solid line in Fig. [T51 shows the results for add- The free-space calculations terminate 
at K = 13.607, where it is known that Efimov physics enters. For k > 13.607, the 
effective dimer-dimer interactions are expected to depend on a three-body or Efimov 
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parameter. 

The dimer-dimer scattering length a^d for unequal-mass systems has also been 
extracted from the energy spectrum of the trapped (2, 2) system with equal frequencies 
and finite-range interactions (see circles in Fig.[T5| |187[|189] . Figure [TBI shows that the 
two sets of calculations (solid line and circles) are in good agreement for n < 13.607. 
The latter calculation does not terminate at k = 13.607 but goes beyond. In particular, 
for K > 13.607, the scattering length add is obtained by following the lowest state of 
the dimer-dimer family, which lies above a subset of states whose energies show a 
pronounced dependence on the three-body parameter. The dimer-dimer scattering 
lengths for k > 13.607 have been confirmed independently by an approximate BO 
treatment (see dotted line in Fig. [16]) [229] . While the scattering lengths derived 
within the BO framework are larger than those obtained by the more exact treatment, 
the overall behavior is similar. A more sophisticated BO treatment, referred to as 
hybrid BO approach, leads to scattering lengths nearly indistinguishable from those 
shown by the solid line and the symbols in Fig. [TBj [229J. The energy spectrum of the 
trapped four-fermion system additionally allows for the determination of the dimer- 
dimer effective range reff,dd (squares in Fig. [T6|) |187[ 1189] . Figure [TH shows that 
''cff.dd increases more rapidly with k than add, indicating that finite-range effects for 
the dimer-dimer system become more important with increasing k. This implies, e.g., 
that the validity regime of the dimer-dimer zero-range model decreases with increasing 
K. The positive dimer-dimer scattering length suggests that two-component Fermi 
gases with large mass ratio k and small atom-atom s-wave scattering length form 
repulsively interacting metastable molecular Bose gases. Losses due to the formation 
of Efimov trimers in dimer-dimer collisions are predicted to decrease exponentially 
with increasing mass ratio |229| . 

Next, we consider the unitary regime and assume that the angular trapping 
frequencies uji are adjusted such that the harmonic oscillator lengths Oho.i, where 
aho.i — ^/fi/inT-iUii), are equal. Small trapped two-component systems with equal 
oscillator lengths, < 20 and 1 < k < 20 have been studied by the fixed node DMC 
method jl211 1189] . These fixed node DMC calculations are performed using a guiding 
function that is constructed from the solution of the two-body problem, a so-called 
"pairing function" . Since the densities of the two components overlap to a fairly large 
degree, the fixed node DMC calculations can be used to estimate the dependence of 
^tr on K. The fixed node DMC energies for k > 1 suggest a slightly modified analysis 
compared to that outlined around Eqs. (|47|) and (|48)) for k = 1. In particular. Eat 
is determined by treating the parameter cetf in Eq- (|48p as a fitting parameter and 
not as a constant |189| . Figure [iTj shows that ^tr (circles in the main figure) and cetf 
(squares in the inset) decrease with increasing k. For comparison, triangles and a 
diamond show the values of ^hom for the homogeneous system obtained by performing 
fixed node DMC calculations [2131 1230] . The slight decrease of ^hom from k = 1 to 6.5 
(triangles) has been interpreted as providing a measure of non- s-wave contributions 
for the K = 6.5 system compared to the k = 1 system [230) . The energy of the 
trapped system (see circles in Fig. [iTl) appears to be lowered more than that of the 
homogeneous system (see triangles in Fig. [T7| . which might be a consequence of the 
fact that the trapped system lowers its energy if the densities are slightly mismatched. 
Other causes for the discrepancy cannot be ruled out and further investigations are 
needed to clarify this point. 

Our discussion thus far has been based on the assumption that trimers are absent. 
From the zero-range model with infinite ao(0), however, it is known that three-body 
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Figure 17. (Color online) Circles (main figure) and squares (inset) show ^tr and 
CbtF) respectively, as a function of the mass ratio k (the values are taken from 
Fig. 6 of Ref. I189| ). For comparison, triangles show ^hom obtained by the fixed 
node DMC method |230| : the diamond shows a more recent value for ^hom for 
K = 1 I2T3I. 



bound states may appear for 8.619 < k < 13.607 [JTJ I161| . The possibility that 
three-body physics, and maybe even A^-body physics, could play a decisive role in 
determining the properties of two-component Fermi gases for sufficiently large n is 
intriguing and has recently stimulated interest among theorists |161[ I163[ I164[ I165| . 
In particular, a three-body bound state with = 1^ symmetry and a four-body 
bound state with = 1+ symmetry have been predicted to appear for the FFX 
systems with k k, 12.31 and the FFFX system with k m 10.4, respectively i 1 Ml HM] . 
These results are obtained by studying trapped systems with purely attractive short- 
range Gaussian model interactions using the stochastic variational approach. An 
independent fixed node DMC study of the free-space system [165j . in which the F 
and X atoms interact through different short-range model potentials, arrives at similar 
conclusions, and furthermore predicts the existence of a resonance for the (4, 1) system 
for K, slightly smaller than 10.4. Although the observed resonances depend, at least 
in principle, on the details of the underlying two-body potential, it has been argued 
that these resonances occur most likely if ~ 1/2 |163[ I164| . The analysis that 
leads to this speculation is based on the hyperspherical framework, which allows for 
the dependence of the eigenenergies on the hyperradial boundary condition to be 
studied. In particular, it was found |163| I164j that the energy- dependence of the 
transcendental eigenequation is weakest for = 1/2 and that the iV-body resonances 
might thus be most amenable to tuning if s^ w 1/2. Although the mass ratio of 
two-component Fermi gases is fixed by nature (implying a particular Si, value), the 
effective mass may be tuned by loading the system into an optical lattice [161] . Future 
work needs to investigate to which extent A^-body resonances can be utilized to study 
or engineer novel many-body phenomena. A first step in this direction was undertaken 
by Nishida et al. |161| who investigated the behavior of many-body systems with two- 
and three-body resonances and k « 8.619 and 13.607. 

Lastly, we consider the regime where n is greater than 13.607, i.e., the regime 
where three-body Efimov physics plays a role (see Subsec. [3?3| . While the three-body 
Efimov effect is quite well understood by now, it is presently an open question to 
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which extent the behavior of the three-body system is modified due to the presence 
of other particles. While such modifications due to the "many-body background" 
are expected to be relatively small for bosons, the situation may be different for 
fermions due to Pauli blocking. In a lowest order semi-classical WKB treatment, 
the presence of the fermions that "do not belong to" the Efimov trimer have been 
accounted for by constraining the momenta hk that contribute to the three-body wave 
function to k > kp [23 Ij . As a consequence, the energy spacing between consecutive 
Efimov states depends, in addition to s,y, on kp- While more studies are needed, 
the outlined example points toward a future frontier, namely the study of few-body 
phenomena in the presence of a many-body background. In passing, we also mention 
that a four-body Efimov effect has been predicted to exist for the FFFX system with 
13.384 < K < 13.607 and 1+ symmetry, i.e., in a regime where the three-body Efimov 
effect is absent |232) . 

4.4- Microscopic approach to thermodynamics of equal-mass Fermi gas 

So far, this review has focused on the zero-temperature regime. An important question 
is, however, how the system properties depend on temperature (see Refs. |223[|233] for 
very recent experiments). Recently, the finite-temperature properties of macroscopic 
two-component Fermi gases have been examined by employing a high temperature 
virial expansion approach that uses the two- and three-body energies as input ^16^ 

[120]. 

For concreteness, we consider a spin- and mass-balanced s-wave interacting two- 
component Fermi gas consisting of N mass m point particles. Since the chemical 
potential /i diverges logarithmically toward — oo as the temperature approaches oo, 
the fugacity z, z = expln/ {kBT)], can be identified as a small parameter |2341l235ll236) . 
It has been shown that a cluster expansion in terms of the fugacity provides reliable 
results for spin- and mass-balanced Fermi gases down to about half the Fermi 
temperature if the expansion includes the second and third order virial coefficients 
(i.e., if the two- and three-body clusters are accounted for) |1691 1170] . The high- 
temperature virial expansion approach is equivalent to an expansion in terms of the 
"diluteness or degeneracy parameter" riA^g, where n denotes the density and AdB 
the de Broglie wave length ^234, ,235. i236i . At high temperatures, the degeneracy 
parameter is small. However, near the transition temperature, i.e., near T/Tc ~ 1, 
the degeneracy parameter is of order one (or AdB ~ (r) , where (r) denotes the average 
interparticle spacing). Thus, the applicability of the virial expansion approach is 
expected to break down in the low-temperature regime. 

Following Liu et al. |169i 1170) , we work in the grand canonical ensemble and write 
the thermodynamic potential in terms of the partition function Z, = —ksTlnZ 
with Z — trexp[— (i7 — fiN)/ (ksT)]. Once the thermodynamic potential il is known, 
other thermodynamic observables such as the energy and pressure can be derived from 
fl by taking appropriate derivatives. The idea behind the virial expansion is to write 
the partition function Z in terms of the partition functions Qn of the n*'' cluster, 
Qn = tr„ exp[— iJ„/(fcBr)] and Z — 1 -\- zQi -\- z^Q\ 4- • • •. Here, H and 7J„ denote 
the Hamiltonian of the full A^-particle system and the n-particle subsystem or cluster, 
respectively. The notation tr„ indicates that the trace is taken over the n-particle 
states. If we denote the eigenenergies of the n*'' cluster by Enj, then the partition 
function Qn becomes Qn — X]j 6^P[~-^n,i/(^s'^)]' where j collectively denotes the 
set of quantum numbers (i.e., the Qn are fully determined by the eigenenergies of 
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Figure 18. (Color online) Shifted virial coefficients Ab„ as a function of the 
dimcnsionless inverse temperature di, di = Hlo / {kBT), of the trapped equal-mass 
two-component Fermi gas at unitarity. Solid and dotted lines show respectively 
the virial coefficient Af)2 and its Taylor expansion up to order Gj^ . Dash-dash- 
dotted and dashed lines show respectively the virial coefficient Aba |237l and its 
Taylor expansion up to order . 



the n*'* cluster). Inserting the expansion for Z into the thermodynamic potential VI 
defines the virial coefficients 6„, Vl = —kBTQi{z + 622:^ + bsz^ + •••). The second 
virial coefficient 62, e.g., equals (Q2 — Qi/2)/(5i and depends on the s-wave scattering 
length ao(0) as well as on the confining geometry and temperature. The expansion of 
in terms of z shows that, if |5„z/6„_i| is smaller than one, the n*'' contribution is 
suppressed by a factor of |6„z/6„_i| compared to the (n — I)*'' contribution. 

We now consider an application to the harmonically trapped two-component 
Fermi system at unitarity, for which the two- and three-particle energies are known 
with high precision |104i If 58| . It is convenient to express the virial coefficients 
6„(ao(0)) of the interacting system relative to the virial coefficients 6„(ao(0) ~ 0) 
of the non-interacting system, i.e., we define A6„ — 6„(ao(0)) — 5„(ao(0) = 0). At 
unitarity, the s-wave two-body energies equal {2j + l/2)hu! + Ecm {j = 0, 17 • • •) 
and A62 can be calculated analytically by summing over the quantum numbers 
associated with the relative and center-of-mass degrees of freedom. The result is 
A62 = exp(— J)/2)/[2(l -I- exp(— J)))] [1691, where Co denotes a dimcnsionless "inverse 
temperature" or scaled frequency, Co = hw / {kBT). The third, shifted virial coefficient 
A63 requires summing the three-body energies, which can be obtained by solving 
a transcendental equation [see Subsec. 13.21 and Eq. (7) of Ref. [158]], numerically. 
The high temperature expansions of A62 and A63 at unitarity are f/4 — w^/32 and 
-0.06833960 + 0.038867J)2, respectively [169]. Figure [H shows the shifted virial 
coefficients A62 and A63 at unitarity as a function of cj. As can be seen, the 
temperature dependence of the virial coefficients is fairly weak in the high temperature 
(small Co) regime. The Fermi energy Ep of the harmonically trapped non-interacting 
two-component Fermi gas, which can be written as Bp ~ {3N)^^^huj, defines the Fermi 
temperature Tp, Tp = Ep/ks- For N = 100, e.g., one finds ksTp w Ihio, i.e., the 
temperature dependence of the virial coefficients is fairly weak if the temperature T 
is larger than the Fermi temperature Tp. 

Knowing the virial coefficients A62 and A&3, the thermodynamic properties of 
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the trapped Fermi gas can be calculated. As a first application, Liu et al. |169| 
determined the interaction energy -Eint, Emt — E — E^J^ (here, E^J^ denotes the energy 
of the corresponding non- interacting system), of a trapped two-component lithium 
mixture at unitarity in the temperature regime T > Tp. The theoretical predictions 
connect smoothly to experimental data from the Duke group |238j . which cover the 
temperature regime T < Tp. More recently, the virial expansion approach has been 
used to determine other observables at and away from unitarity [239. .240 . 

The virial expansion approach outlined can be extended to the homogeneous 
system, thereby allowing for comparisons with results obtained by finite temperature 
MC simulations |2411 1242i 1243] and effective many-body approaches. To this end, Liu 
et al. |169j applied an "inverse LDA" , which allowed for the determination of the virial 
coefficients of the homogeneous system from those of the inhomogeneous system. The 
theoretically determined high-temperature value of the third order virial coefficient 
has been confirmed experimentally [223) . The examples discussed in this subsection 
underline that microscopic bottom-up approaches provide, in certain cases, powerful 
means to predict many-body properties. 

5. s-wave interacting Bose gas under confinement 

This section considers single-component Bose gases under external confinement. 
Subsection 15.11 introduces the mean-field Gross-Pitaevskii (GP) framework, which 
assumes that the behavior of weakly-interacting Bose gases is fully governed by the s- 
wave scattering length ao(0). The validity regimes of the GP equation and a modified 
GP equation are benchmarked for Bose gases under spherically symmetric harmonic 
confinement. Subsection 15.21 discusses the mechanical instability that arises for Bose 
gases with negative s-wave scattering lengths. Lastly, Subsec. lOl discusses the role of 
the effective range and the three-body parameter for Bose systems. 

5.1. Equation of state of trapped Bose gas: The role of the s-wave scattering length 

Weakly-interacting Bose gases under external confinement at zero temperature are 
most commonly described by the non-linear mean-field GP equation [TTJ [TH [131 111 
I244[ 12451 1246j , which can be derived by approximating the many-body wave function 
\l/(ri,- • • ,f]\f) by a product of single-particle orbitals ^{rj), 4* — 11^=1 ^{^j)^ where 
/|$(r)pd^r = 1 [ini 1247) . This product ansatz is fully symmetrized and geared 
toward the description of the energetically lowest lying gas-like state of the iV-boson 
system. Assuming that the interactions can be parameterized through a sum of two- 
body potentials Vthir jk) and minimizing the total energy Etot{N) with respect to $*, 
one obtains the Hartree equation [TH 1247) 



In Eq. (|50p. m denotes the atom mass, Vtrap the external confining potential 
felt by each of the bosons, and e the orbital energy or chemical potential. The 
integro-differential equation (I5(I1) . which has to be solved self-consistently for 4>(r) 
and e, describes a trapped "test particle" that moves — ^just as in atomic structure 
calculations — in the mean-field [the third term in square brackets on the left hand 
side of Eq. ([50]) ] created by the other {N — 1) atoms. 




e$(r). 



(50) 
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For ultracold degenerate s-wave interacting Bose gases, the two-body interaction 
potential Vtb with s-wave scattering length ao(0) is frequently replaced by the Fermi 
pseudopotential Vp, Eq. for which the integral in Eq. ((5(7)) can be evaluated 

analytically, 



2^-V^^+ FtrapW -f (TV - 1) 



e$(f). (51) 

Equation (|5ip is the GP equation for s-wave interacting systems, sometimes also 
referred to as non-linear Schrodinger equation. Although the GP equation has been 
derived by applying the variational principle, the resulting total energy Etot{N) does 
not constitute a variational upper bound to the energy of the energetically lowest- 
lying gas-like state of the many-body system since the true atom-atom potential has 
been replaced by Fermi's pseudopotential. The GP equation depends on the product 
{N — l)ao(O), and not on (TV — 1) and ao(0) separately. This implies that the mean- 
field formulation is only sensitive to the "strength" of the effective interaction and not 
to its underlying microscopic origin, i.e., whether it is due to, say, large N and small 
ao(0) or small N and large ao(0). An alternative derivation of Eq. (ISTI) is based on a 
field operator description and results in a factor of TV as opposed to TV — 1 (see, e.g., 
Refs. [H mg). WhUe it is justified to replace TV - 1 by TV in the large TV limit, 
the factor of TV — 1 should be used in the small TV limit where the difference between 
TV — 1 and TV is appreciable. The outlined number conserving microscopic derivation 
of the GP equation emphasizes that its application is not restricted to large number of 
particles; in fact, as shown below, the GP equation provides a good description even 
for dilute two-particle systems. 

The sohd line in Fig. [TOl shows the total energy Etot/N per particle, obtained 
by solving the GP equation for a spherically symmetric harmonic confinement, as a 
function of the interaction parameter (TV — l)ao(0)/aho- The solid line terminates at 
(TV — l)ao(0)/aho ~ —0.57497 |247[ I250j . For more negative interaction parameters, 
the GP equation supports a solution whose energy is unbounded from below (i.e., 
the energy approaches —oo and the corresponding wave function describes a high 
density state of size much smaller than Oho)- Generally speaking, the validity regime 
of the GP equation is restricted to dilute systems, i.e., to parameter combinations that 
fulfill the inequality n(0)|ao(0)|'^ <C 1 |251| . where n(0) denotes the peak density. For 
(TV — l)ao(0)/aho ^ —0.57497, the diluteness parameter n(0)|ao(0)P is much smaller 
than 1, assuming TV is not too small. This suggests that the predicted drop of the 
GP energy (the termination of the solid line in Fig. [T5| lies within the validity regime 
of the mean-field framework. As discussed in Subsec. 15.21 the sudden drop of the GP 
energy is interpreted as a mechanical instability or collapse. 

We now compare the solutions of the GP equation with those obtained by solving 
the linear Schrodinger equation. Dotted, dashed and dash-dotted lines in Fig. [T9l show 
the total energy Etot/N per particle for two particles under spherically symmetric 
harmonic confinement interacting through a hardsphere potential with range ao(0) 
[ao(0) > 0], a square well potential with range ro and depth Vq (ro = O.Oloho and Vq 
are chosen such that the potential supports one s-wave bound state in free space), and 
the pseudopotential Vq^ [see Eq. (jj])], respectively. The two-body energies for these 
potentials are obtained by solving the linear Schrodinger equation analytically (for the 
treatment of Vq^, see Sec. [5]). For (TV — l)|ao(0)|/aho ^ O.f , the agreement between 
the energies obtained by solving the linear and non-linear Schrodinger equations is 
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Figure 19. (Color online) Total energy Etot/N per particle of a dilute Bose gas 
under external spherically symmetric confinement as a function of the "interaction 
parameter" (N — l)ao(0)/aiio. The solid line shows the mean-field GP energy. 
Dotted, dashed and dash-dotted lines show the energies obtained by solving the 
linear Schrodinger equation for two particles interacting through a hardsphere 
potential, a square-well potential (with range ro, ro = O.Olaho, and depth Vq 
adjusted such that the free-space system does support one bound state) and the 
zero-range pseudopotential Vg''' [Eq. l(4]l], respectively. Note that the dashed and 
dash-dotted lines arc indistinguishable on the scale shown. 



excellent. However, for larger (A^— l)|ao(0)|/aho deviations are visible. The energy for 
two particles interacting through the hardsphere potential increases roughly linearly 
with ao(0) and lies above the other energies. The comparatively large increase of 
the energy for the hardsphere potential can be understood by realizing that the 
excluded volume increases with increasing ao(0), resulting in an infinite two-body 
energy as ao(0) approaches infinity. This unphysical behavior suggests that the 
applicability regime of the hardsphere potential is limited to small ao(0). The GP 
energy lies below the energies for the zero-range and square- well interaction potentials 
for {N — l)ao(0)/aho ^ 1 and above for {N — l)ao(0)/aho 1. The discrepancy 
between the GP energy and the energies for the zero-range and square-well interaction 
potentials increases as {N — l)ao(0)/aho increases further. In the large interaction 
parameter regime, the diluteness parameter n(0)|ao(0)p is no longer small compared 
to one and higher order corrections need to be included in the mean-field description 
(see below). For negative scattering lengths, the energies obtained by solving the 
linear Schrodinger equation vary smoothly. In particular, the mechanical instability 
predicted by the GP equation has no two-body analog; at least three particles are 
needed to see an analog of the instability predicted by the GP equation within the 
framework of the linear Schrodinger equation (see also Subsec. [221). 

For larger number of particles, the determination of the solutions to the linear 
Schrodinger equation becomes more involved. In addition to the increased number 
of degrees of freedom, complications arise due to the fact that the iV-boson system 
supports, unlike s-wave interacting two-component Fermi systems, A^-body (N > 2) 
bound states even if the underlying two-body potential does not support a two- 
body bound state. Thus, the energetically lowest lying gas-like state of the Bose 
gas is a metastable state and not the true ground state of the system. Tightly 
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bound A^-body clusters, whose properties are discussed further in Subsec. 15. 3[ are 
absent for systems interacting through purely repulsive two-body potentials such as 
the hardsphere potential. Although purely repulsive interaction models exhibit, as 
discussed exemplarily above for the hardsphere potential, unphysical behavior for 
large ao(0), they provide a quantitatively correct description for small oo(0) and allow 
one to estimate model dependent effects for intermediate ao(0) |2521 12531 1254[ [255] . 
Symbols in Fig. [20l a) show the energy Etot/N per particle for A^-boson systems, N = 
2 — 50, with hardsphere interactions under external spherically symmetric harmonic 
confinement as a function of the interaction parameter {N — l)ao(0)/aho [252) . It 
can be seen clearly that the solutions of the linear Schrodinger equation depend on 
N and ao(0)/aho separately and not just the product of these two parameters. For 
comparison, the solid line shows the GP energy per particle. The DMC energies for 
different number of particles (symbols), obtained by solving the linear Schrodinger 
by the DMC method [252j . lie above the GP energies. Figure [^CTb) shows the energy 
difference E'dmc— ^'GP- For fixed N, the deviations between the DMC and GP energies 
increase with increasing interaction parameter. For a fixed interaction parameter, the 
deviations decrease with increasing N. An improved mean-field description of the 
hardsphere equation of state is discussed below. 

To investigate the small |ao(0)|/aho behavior in more detail, the iV-boson 
system under external spherically symmetric confinement with two-body zero-range 
interactions can be treated perturbatively within the framework of linear Schrodinger 
quantum mechanics. Up to order [ao(0)/aho]^, the perturbative energy of the 
energetically lowest-lying gas-like state of the A^-boson system reads |256j 
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The first term on the right hand side is the energy of the non-interacting (unperturbed) 
A^-boson system. The leading order energy shift of the A^-particle system is equal to 
the energy shift that A^(A^ — l)/2 independent two-body systems would experience. 
The next order shift contains a positive term that scales with the number of pairs 
and a negative term that scales with the number of trimers. The latter contribution 
has been interpreted as an effective attractive three-body interaction |1601 1256) . This 
interpretation emerges naturally in an effective field theory treatment [256) . which 
derives an effective low-energy Hamiltonian that contains one-body, two-body and 
three-body terms. We stress that the effective attractive three-body interaction 
arises from the pairwise zero-range interactions and not from an explicit three-body 
potential. The need for and role of an actual three-body potential is discussed 
in Subsec. 15.31 The effective three-body interaction for the harmonically trapped 
Bose gas was first predicted by Jonsell et al. [160) , who employed the BO and 
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Figure 20. (Color online) Hardsphere equation of state for the Bose gas under 
spherically symmetric confinement. Symbols in panel (a) show the total energy 
EoMc/^ per particle for the A^-particle Bose gas with hardsphere interactions 
as a function of the interaction parameter {N — l)ao(0)/aiio for N = 2 (circles), 
= 3 (squares), N = 5 (diamonds), Af = 10 (up triangles), Af = 20 (left 
triangles), and Af = 50 (down triangles); the energies are obtained by solving 
the linear Schrodinger equation by the DMC method. For comparison, the 
solid line shows the total energy per particle Eqp/N obtained by solving the 
mean-field GP equation. Panels (b) and (c) show the scaled energy differences 
{Et)mc-Egp)/Ey)mc and (-Edmc - -EGP,mod)/-EDMC, respectively, for the same 
A'^ as panel (a); to guide the eye, dotted lines connect data points for the same 
A^. The figure has been adapted from Ref. [252^. 
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adiabatic approximations within the hyperspherical framework to place bounds on 
the three-body coefRcient. For a Bose gas confined to a box with periodic boundary 
conditions, the effective three-body interaction has a positive coefficient ^T] (see also 
Refs. |257[ I258T 12591 f25D] ) . which illustrates that the effective A^-body interactions can 
be tuned by varying the external confinement. The second term on the right hand 
side of Eq. ()52|) can also be derived within the GP framework by approximating $ in 
Eq. ()5ip by the non-interacting wave function. 

Equation (|52p can be viewed as the few-particle trap analog of the low-density 
expansion of the equation of state of the homogeneous Bose gas with density 



Application of the LDA to the first term on the right hand side of Eq. ([55]) gives 
the second term on the right hand side of Eq. (1521) . The next order correction in 
Eq. ([5^ . in contrast, cannot be determined by applying the LDA to the na^{Q) 
term in Eq. ([55]) [160]. 

Lastly, we comment on the large scattering length regime. The experimental 
realization of cold atomic gases with large scattering length and sufficiently long 
lifetimes would open the possibility to study strongly-correlated Bose systems with 
unprecedented control. Historically [2621 1263) , the study of liquid helium, which is 
characterized by a gas parameter naQ(O) of about 0.21 and a condensate fraction 
of about 7% [2641 I265i 1266] . is tremendously important. Liquid helium has an 
s-wave scattering length that is roughly 10 times larger than the effective range, 
which in turn is of the order of the interparticle spacing. Although the s-wave 
scattering length of Bose gases can, for a subset of species such as ®^Rb |33], be 
tuned over a wide range through the application of an external magnetic field in 
the vicinity of a Fano-Feshbach resonance, Bose gases with large scattering length 
exhibit detremental losses since the three-body recombination rate shows an overall 
|ao(0)|'* scaling [TTT I [1721 [1731 [T74 l [267] . This is in contrast to two-component Fermi 
gases with equal masses, which are stabilized by the "Pauli pressure" that results 
from the Pauli exclusion principle. While it is not clear at present to which extent 
the equilibrium properties of Bose gases with large scattering length can be probed 
experimentally [268j . several theoretical studies have investigated the properties of 
Bose gases with large scattering length based on the assumption that the system's 
lifetime would be sufficiently long to reach equilibrium. The first studies in this 
direction were performed for the trapped three-boson system by Jonsell et al. |160| 
and Blume et al. [166] . More recent studies have estimated the lifetime of the 
trapped three-boson system [1 581 1179) . Extensions to larger systems have also been 
considered [2181 1269, 270, 271] |272l [273]. The work by Song and Zhou 272 , e.g., 
predicts the existence of a fermionized three-dimensional Bose gas with large scattering 
length (see also Ref. )274) ). This is an intriguing prospect since the fermionization of 
Bose gases has been studied extensively in one-dimensional systems. 

If the interaction parameter {N ~ l)ao(0)/aho becomes appreciable but not too 
large, corrections to the GP equation can be accounted for by a modified mean-field 



n [2581 [26T] . 




(55) 



GP equation 2521 l275l [2761 [277] . 
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[i + ^[aimN - = ^w)- (56) 



The second term in the big round brackets in the second hne of Eq. (j56|) can be derived 
from the "quantum fluctuation term" of the homogeneous system, i.e., the second 
term on the right hand side of Eq. ([55]) . Figure ISOT c) shows the energy difference 
EjyMC ~ -EGP^mod, where the energy i?GP,mod has been derived from the orbital energy 
e that resuhs when solving Eq. ([55)) . A comparison of the vertical scales of Figs. l^Uf b) 
and [20lc) shows that the modified GP equation extends the validity regime of the 
mean-field treatment. 



5.2. Bose gas with negative s-wave scattering length 

To shed light on the mechanical instability for Bose gases with negative scattering 
length (i.e., the termination of the solid line in Fig. for negative scattering 
length), we solve the GP equation for a spherically symmetric harmonic confinement 
using a variational Gaussian wave function with width b. Writing $(r^ = 
6~^/^7r~^/'* exp[— r^/(26^)], the variational energy becomes |278j 



N 



3aL _^ 36^ , {N-l)ao{0)at 



huj, (57) 



462 4a2„ V^fe3 

where the first, second and third terms on the right hand side are the kinetic, 
trap and interaction energies, respectively. Since the interaction energy varies as 
b~^ while the kinetic and potential energies vary as b~^ and b^, respectively, it 
is clear that there exists a critical negative interaction parameter beyond which 
the variational energy possesses a global minimum for b = but not a local 
minimum |278| . Similar conclusions have been drawn in related studies [2791 12801 1281j . 
Lines in Fig. I21f a) show the variational energy per particle as a function of 6/aho 
for various interaction parameters {N — l)|ao(0)|/aho- For small {N — l)|ao(0)|/aho, 
the variational energy shows a minimum around b/uho ~ 1- As (N — l)|ao(0)|/aho 
increases [ao(0) negative], the energy minimum moves to smaller b/uho and the 
barrier that separates the local and global energy minima decreases. At the critical 
value of {N — l)ao(0)/aiio ~ —0.67051 [see dashed line in Fig. [2]Ta)]. the barrier 
disappears, suggesting that the Bose gas can lower its energy — without having to 
"tunnel" through an "energy barrier" — by shrinking toward (or collapsing into) a so- 
called high density "snowfiake state". The variational Gaussian approach describes 
the onset of the mechanical instability quite well, i.e., it predicts a slightly more 
negative critical interaction parameter than the full solution to the GP equation 
[{N - l)ao(0)/aho ~ -0.67051 compared to -0.57497]. Furthermore, it highlights 
that the trapped Bose gas is, for sufficiently small (TV — l)]ao(0)l/aho [ao(0) < 0], 
stabilized by the positive energy contributions that arise due to the presence of the 
trap. If the system was homogeneous, an infinitesimally small |ao(0)] [ao(0) < 0] 
would induce collapse. While it is suggestive to interpret the energy curves shown in 
Fig.mja) as effective potentials, it must be kept in mind that the width b corresponds 
to a variational parameter and not to an actual spatial coordinate. 

The collapse dynamics of Bose gases with negative oq (0) has also been investigated 
by solving the linear Schrodinger equation with zero-range interactions within the 
BO hyperspherical approximation [149[ . For a variational wave function that is 
written as the product of a constant hyperangular piece and a hyperradial function 
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R/a, 



ho 



Figure 21. (Color online) Mechanical instability of Bose gases with negative 
scattering length, (a) Dotted, dash-dot-dotted, solid, dashed and dash-dotted 
lines show the variational energy Evar/N per particle, Eq. II57I I. as a function 
of b/a^^a for (A'^ — l)ao(0)/aijo = 1,0,-0.3,-0.67 and —1, respectively, (b) 
Dotted, dash-dot-dotted, solid, dashed and dash-dotted lines show the effective 
potential curve Vcg{R) / N , Eq. jsHJl, as a function of R/a-^^ for {N—l)aQ (0)/aho = 
1, 0, -0.3, -0.67 and -1, respectively [N is fixed, N = 1000). 



R ^^/'^F{R), the effective fiyperradial potential curve Vcs{R) takes tlie form |149] 

n^{3N -1){3N -3) 1^^ 

f ) = SnIr^ + 2^"^- ^ + 

1 n^(iV-l)ao(0) r(3jV/2) N 

2tt m r((37V-3)/2)iV3/2i?3' ^ ) 

wliere NR^ = The effective potential V;,ff(i?) enters into the hyperradial 

equation —h^/{2Nm)F" + VcgF = EtotF, which can be solved readily. In Eq. (|58|) . 
the first, second and third terms on the right hand side arise due to the kinetic energy 
operator, the external confinement, and the pairwise interactions, respectively. Lines 
in Fig.lCTb) show Vce{R)/N for various interaction parameters {N — l)ao(0)/aho and 
N = 1000. We note that the curves Vcs{R)/N depend only weakly on TV [T^ . A 
comparison of Figs. dHa) and YlVih) [and, equivalently, of Eqs. ((57)) and ([55]) ] shows 
that b plays a role similar to that played by R. The key difference is that 6 is a 
variational parameter while the hyperradius i? is a spatial coordinate. Physically, the 
small R region describes a system in which the atom-atom distances are small, i.e., 
where the system is self-bound. Using the effective potential curves T4ff, tunneling 
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from the region where the potential exhibits a local minimum to the small R region 
where the effective potential is unbounded from below has been estimated within the 
WKB approximation 11491 . The conclusions discussed here have been confirmed for 
the three-body system by performing an analysis that goes beyond the variational 
treatment and furthermore accounts for the coupling between channels [166) . For 
larger systems, confirmation comes from treatments that explicitly account for two- 
body correlations |282] and experiments [283, 284,. 

5.3. Beyond the s-wave scattering length: Effective range and three-body parameter 

The previous two subsections assumed that the equation of state of trapped Bose 
gases is fully determined by the s-wave scattering length ao(0). However, atom-atom 
potentials possess a finite range, which induces range-dependent corrections to the 
equation of state. Furthermore, the properties of Bose gases may additionally depend 
on a three-body parameter whose role for the three-boson system with large s-wave 
scattering length has already been discussed in Subsec. 13.31 In the following, we 
present simple arguments that illustrate at which order these corrections appear in 
the equation of state. 

To account for the finite range of the underlying atom-atom potential, we resort 
to the effective range expansion, Eq. The effective range rcs is determined by 
the shape and range of the two-body potential. For the hardsphere potential, e.g., 
the effective range rcfi is directly proportional to ao(0) [i.e., Tcff — 2ao(0)/3]. For the 
van der Waals potential with r~^ tail, the effective range r^g is, in the low-energy 
regime, also determined by the scattering length ao(0) [285]. To estimate the order 
at which finite range effects occur, we replace the s-wave scattering length ao(0) in 
the leading order term of Eq. (j52p by ao(-E) and, in a second step, approximate ao{E) 
by Eq. ([3]). A simple calculation shows that finite range corrections enter, according 
to this estimate, at order rcffaQ(0)/aj^Q. Thus, finite range effects are, assuming that 
I Toff I is smaller than Oho, suppressed compared to the effective three-body interaction 
[see Eq. (|52|) in Subsec. 15. 1| . which scales as [ao(0)/aho]^, by a factor of rctf/oho- 
The calculation could be made rigorous by carrying out a perturbative analysis for 
a sum of zero-range two-body potentials that account for finite-range effects through 
derivative operators 81 [i.e., that depend on ao(0) and r^s]. Within the mean-field 
framework, effective range corrections have been accounted for by adding a gradient 
term that is proportional to g2{N — l)V^|$(r)|^ to the mean-field potential of the 
GP equation 286 . To determine this modified mean- field term and the constant g2, 
Fu et al. derived a pseudopotential that reproduces the real part of the scattering 
amplitude up to order 12861 and then used this revised pseudopotential in the 
Hartree framework (see Subsec. [O]) . 

The role of the three-body parameter has previously been discussed in Subsec. [3.31 
for the three-boson system with infinitely large two-body scattering length ao(0). We 
now consider the small ao(0) regime and estimate at which order the three-body 
parameter enters into the equation of state. As already discussed, an imaginary Si, 
gives rise for the need of a three-body parameter. For small |ao(0)|. Fig. [5] shows that 
there exists an imaginary root for sufficiently small R. This implies that the sum of 
zero-range two-body potentials has to be complemented by a hyperradial three-body 
potential V3 or, equivalently, a three-body boundary condition in the hyperradial 
coordinate. While this boundary condition does not, as in the large |ao(0)| limit, lead 
to the appearance of an infinite ladder of geometrically spaced finite energy states, it 
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does have an effect on the energy spectrum. 

To estimate the effect of V3 on the energy spectrum of the trapped three-boson 
system, we perform a perturbative analysis. We assume that V3 is proportional 
to {h?A'^{k)/m)6{ri — r2)5{r2 — rs), where A'^{k) denotes the generalized energy- 
dependent three-body scattering length. Dimensional analysis shows that A'^{k) has 
units of (length)'*. This three-body pseudopotential can alternatively be written 
in terms of a 5-function in the hyperradius. The quantity ^^(0) conveniently 
parameterizes the three-body parameter in the weakly-interacting regime and is 
related to the hypervolume D defined by Tan f260| . In first order perturbation theory, 
this pseudopotential leads to an energy shift proportional to [yl(fc)/aho]^fiw. For the 
A'^-boson system, this energy correction gets multiplied by the number of trimers. This 
analysis reveals that the three-body potential introduces an energy shift of higher order 
in the inverse harmonic oscillator length than both the effective three-body interaction 
and the two-body effective range correction. We note that the effect of the three-body 
parameter on the energy spectrum of the weakly-interacting iV-boson system in a box 
with periodic boundary conditions has been discussed in Refs. |259i 1260] . 

We now turn to the large ao(0) regime and ask two questions: (i) Does the 
description of successively larger Bose systems require a new parameter for each 
successively larger weakly-bound Bose system? (ii) What are the imprints of the three- 
body parameter on the A^-body {N > 3) spectrum? The first question has been studied 
extensively ever since Efimov's work from the early seventies [45l |46[ l47l I287i [2881 1289j . 
It is fairly well established by now that the treatment of the low-energy physics of 
A^-body {N > 3) Bose systems does not require additional parameters, i.e., the low- 
energy properties of the iV-body system are believed to be determined by ao(0) and 
a three-body parameter |289[ 12911 [292j . A consequence of the absence of four- and 
higher-body parameters is that there exists no four- or higher-body Efimov effect 
for Bose systems with 0"*" symmetry, i.e., there exists no infinite sequence of iV-body 
bound states whose geometric spacing is determined by an A^-body parameter (A^ > 3). 
Other interpretations do, however, exist [290]. The second question has been studied 
extensively over the past five years or so, mostly but not exclusively in the context of 
the four-body system in free space |2911 1292] . For each Efimov trimer, there exist two 
universal four-body states whose properties are determined by those of the Efimov 
trimer that the four-body states are "attached to". These four-body states are, in 
fact, resonance states with finite lifetime |295j . The energies of the tetramer "ground 
state" and "first excited state" in free space are about 4.58 and 1.01 times larger 
than the binding energy of the respective trimer |2911 1292) . The factor of 4.58 can 
be rationalized by realizing that the tetramer can be thought of as consisting of four 
trimers. In addition to the energy spectrum, other observables — such as the scattering 
lengths for which the tetramers become unbound and dissociate into four free atoms — 
obey universal relations |293| I294[ I295[ I296| [54 ] . At present, the study of larger systems 
is still in its infancy. Initial predictions that relate the A^-body energies to those of 
the trimer have been made |289[ 12971 12981 1299] and it is expected that this research 
area will florish in the years to come. 

6. Summary and Outlook 

Trapped ultracold atomic and molecular gases present a research area that is 
blossoming and whose future looks bright. This area has attracted the attention 
of researchers from diverse areas, including those with backgrounds in atomic physics. 
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nuclear physics, few-body physics, quantum optics, condensed matter physics and high 
energy physics. One attractive aspect of this research area is that experimental and 
theoretical efforts go hand in hand, with experiments being motivated by theoretical 
work and theoretical work being motivated by experimental progress. 

This review provides an introduction to trapped atomic and molecular physics 
following a bottom-up approach. Building on a detailed discussion of two-body and 
three-body systems, the properties of trapped fermionic and bosonic systems with 
varying number of particles have been investigated. The discussions and examples 
presented are meant to give the reader a flavor of this rich and rapidly developing 
interdisciplinary field. In the following, we comment on selected present and future 
frontiers that have only been touched upon or not been covered at all in this review. 

Other mixtures: This review focused primarily on the properties of two- 
component Fermi gases and single-component Bose gases. Extensions of these studies 
to other mixtures, such as Bose-Fermi or Bose-Bose mixtures, introduce additional 
degrees of freedom (i.e., another scattering length and possibly another confining 
length and/or mass ratio) and will allow one to investigate the interplay between 
interactions and symmetry in more depth. Studies along these lines may, e.g., lead to 
a bottom- up understanding of induced interactions (see, e.g., Ref. |300j ). 

Multi- component Fermi gases: While two-component Fermi gases are stable for 
essentially any interaction strength, Fermi gases with more than two components are 
not necessarily stable against collapse. In fact, the phase diagram of multi-component 
Fermi gases is currently being debated |30H 13021 i303, 304] I305j . Bosenova-type 
collapse, BCS-type phases and phases whose properties are governed by competing 
dimer and trimer formation have been suggested. Studies of trapped multi-component 
Fermi gases may shed light on the effective interactions between two trimers and its 
connections to Efimov physics. 

p-wave Fermi gases: This review focused almost exclusively on s-wave interacting 
systems in the vicinity of broad Fano-Feshbach resonances. While present experiments 
on higher partial wave resonances are plagued by detrimental losses, it seems 
feasable that these losses could be suppressed by loading the gas into an optical 
lattice or by utilizing atomic species for which dipolar relaxation rates can be 
suppressed [306] . Some theoretical few-body studies on p-wave interacting systems 
have been performed j307| I308| 13091 1310j , and it is anticipated that additional studies 
will further clarify questions related to the system's lifetime and address which aspects 
of higher partial wave systems are universal and which ones are not. 

Dipolar systems: Much progress has been made in trapping and cooling systems 
with magnetic and electric dipoles. While weakly-interacting dipolar Bose gases have 
been modeled quite successfully using a mean-field description, the study of strongly- 
interacting dipolar systems is still in its infancy. The three-body system has been 
investigated for a few limiting scenarios within a microscopic framework |31H I312j 
but little has been done yet for larger systems. One possible approach is to pursue 
a description that is based on an effective Hamiltonian, which incorporates our 
understanding gained from microscopic two- and three-body studies. 

Extended Efimov scenario and N-body parameter: The Efimov effect, in its 
original formulation, applies to three-body systems in free-space. Ever since Efimov's 
papers from the early 1970s, researchers have asked whether or not the Efimov effect 
generalizes to larger systems. Although much progress has been made, there are 
still open questions. Future studies of free-space and trapped four-, five- and higher- 
body systems will address in more detail how the energy spectrum of the iV-body 
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system [N > 3) depends on the three-body parameter and how finite-range corrections 
enter. It is anticipated that answers to these questions will be found by combining 
the insights gained from microscopic and effective approaches applied to systems of 
varying symmetry and size. 

Systems with two- and/or three-body resonances: This review discussed how 
two- and three-body interactions can be tuned by varying the trappping frequency 
of the harmonic confinement, by taking advantage of the optical lattice structure, 
or by inducing tunable eflFective three-body interactions that correspond to varying 
short-range hyperradial boundary conditions. For example, the scattering properties 
of three- and four-body systems in effectively one- and two-dimensional confining 
geometries have been analyzed in some detail |313[ 1314] . Another intriguing example 
is the appearance of resonances in mixed dimensional systems |315[I316] . It is expected 
that future work will uncover other means to manipulate the interactions. The 
resulting effective interactions are expected to open unique opportunities for designing 
larger systems with exotic properties. This is a research branch that appears to be 
equally intriguing to few- and many-body physicists. 

Input for Hubbard Hamiltonian: Any type of effective Hamiltonian such as the 
Bose or Fermi Hubbard Hamiltonian depends on effective parameters such as the 
tunneling amplitude and on-site interactions. As experiments are being refined and 
many-body simulation techniques advance, it is paramount to refine the determination 
of the effective parameters that define a given effective Hamiltonian |3171I318] . Studies 
of trapped few-body systems are expected to contribute to this endeaver. 

Beyond the LDA: Connections between inhomogenous and homogeneous systems 
are often made through the LDA. While this approach has been tremendously 
successful, several limitations of this approach are known. A future goal of trapped 
atom and molecule studies should be to develop alternative approaches. These could 
involve pushing microscopic calculations of trapped atom systems to larger number of 
particles and extending condensed matter techniques to trapped systems. 

Time- dependent dynamics of Bose gases: Although this review focused primarily 
on time-independent studies, it should be noted that much progress has been 
made over the past few years in solving the time-dependent many-body Schrodinger 
equation. One-dimensional Bose gases, e.g., have been treated by various variants 
of the multi-configurational approach |319) (borrowed from quantum chemistry) 
and the time-evolving block decimation or density-matrix renormalization group 
approaches |320| I32H 1322] (borrowed from condensed matter physics). From a 
fundamental as well as an experimental point of view, the emergence of fragmented 
states of Bose gases is particularly interesting. Another intriguing set of studies 
addresses the tunneling dynamics in double-well geometries. In particular, the 
interplay between single particle and pair tunneling has been investigated [323] . It is 
expected that studies along these lines will be extended to higher-dimensional systems 
in the future. 

Finite-temperature physics: Except for Subsec. 14. 4i this review considered system 
properties at zero-temperature. Many features, such as the shell structure of two- 
component Fermi gases [324] and peaks in the loss rate due to Efimov resonances j325j , 
become smeared out as the temperature increases. At the microscopic level, few 
studies to date have considered thermally averaged quantities, despite the fact that 
thermometry in optical lattice systems is currently a hot topic [326] . It is expected 
that the few-body community will push to develop the machinery that allows for the 
treatment of trapped systems at finite-temperature. 
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N on- equilibrium few-body dynamics: One of the key advantages of small trapped 
systems is that their entire excitation spectrum is, in certain cases, amenable to 
analytical or numerical treatments. This opens the possibility to prepare a system 
in a given initial state and to then follow the dynamics in response to a slow variation 
of a system parameter or a rapid quench. Initial studies along these lines have been 
performed [72j I188j . It is expected that this area of research will further intensify 
as studies of the non-equilibrium dynamics of optical lattice systems are progressing 
dramatically. Answers related to equilibration times, entanglement dynamics, and 
critical behavior are expected to emerge. For example, it is anticipated that few-body 
studies will address what happens when two initially separated optical lattice sites or 
plaquettes are being merged. 
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